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The  localized  orbitals  discussed  here  will  be  divided  into  two 
classes:  (1)  intrinsically  localized  orbitals,  where  the  localization  is 
due  primarily  to  symmetry  or  energy  considerations,  for  example 
transition  metal  d-orbitals  or  lanthanide  /-orbitals;  and  (2)  orbitals 
which  must  be  localized  after  a  self-consistent  field  (SCF)  calculation. 
In  the  latter  case,  two  new  methods  of  localization,  the  Fermi  and  the 
double  projector  methods,  are  presented  here.   The  Fermi  method  provides 
a  means  for  the  non-iterative  localization  of  SCF  orbitals.  while  the 
double  projector  allows  one  to  describe  what  atomic  functions  the 
localized  orbitals  will  contain.   The  third  localization  procedure 
described  is  the  second  order  Boys  method  of  Leonard  and  Luken.   This 
method  is  used  to  explain  the  photodissociation  products  of  Ni(CO).. 
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The  Intermediate  Neglect  of  Differential  Overlap  (INDO)  method  is 
extended  to  the  /-orbitals,  and  the  intrinsic  localization  of  the  /- 
orbitals  is  examined.   This  extension  is  characterized  by  a  basis  set 
obtained  from  relativistic  Dirac-Fock  atomic  calculations,  and  the 
inclusion  of  all  one-center  two-electron  integrals.   Applications  of 

this  method  to  the  lanthanide  halides  and  the  twelve  coordinate 

_2 
Ce(N02)g  ion  are  presented.   The  model  is  also  used  to  calculate  the 

ionization  potentials  for  the  above  compounds.   Due  to  the  localized 

nature  of  /-orbitals  the  crystal  field  splittings  in  these  compounds  are 

extremely  small,  leading  to  SCF  convergence  problems  which  are  addressed 

here.   Even  when  the  SCF  has  converged,  a  small  configuration 

interaction  (CI)  calculation  must  be  done  to  insure  that  the  converged 

state  is  indeed  the  lowest  energy  state.  The  localized  nature  of  the  f- 

orbitals  in  conjunction  with  the  double  projector  localization  method 

may  be  used  to  isolate  the  /-orbitals  in  order  to  calculate  only  a  CI 

restricted  within  the  /-manifold. 
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INTRODUCTION 


Localized  orbitals  may  be  defined  as  either  orbitals  which  are 
spatially  compact  or  as  molecular  orbitals  vhich  are  dominated  by  a 
single  atomic  orbital.  The  use  of  the  terms  bond,  anti-bond,  or  lone 
pair  to  describe  a  set  of  orbitals  are  all  based  on  a  localized  orbital 
framework.  The  use  of  ball  and  stick  models  and  hybrid  orbitals  in  every 
general  chemistry  class  illustrates  the  power  of  localized  orbitals  as 
an  aid  in  the  understanding  of  molecular  structure. 

Localized  orbitals  may  be  divided  into  two  categories.  The  first 
category  encompasses  orbitals  which  must  be  localized.  Although,  some 
orbitals  are  localized  automatically  either  by  their  symmetry  or  by 
their  energy  relation  to  other  orbitals  in  the  molecule,  these  orbitals 
form  the  basis  for  the  second  category  of  localized  orbitals. 
Transition  metal  d-orbitals  fall  into  this  second  category,  and  it  has 
been  predicted  that  the  lanthanide  /-orbitals  should  also  fall  into  this 
category. 

Our  understanding  of  transition  metal  chemistry  is  also  based  on  the 
concept  of  orbitals  being  localized.   The  excitations  that  give  rise  to 
the  colors  of  many  metal  complexes  are  classified  as  d-d,  ligand-d  or 
charge  transfer.   These  classifications  are  based  on  the  fact  the  d 

orbitals  are  localized  allowing  for  the  easy  interpretation  of  spectra. 

The  success  of  crystal  field  theory  reinforces  the  belief  that  the  d- 

orbitals  are  localized. 

In  the  following  two  chapters,  I  will  examine  both  of  these  types  of 

localized  orbitals.   The  first  chapter  will  deal  with  methods  developed 

to  obtain  localized  orbitals  from  delocalized  orbitals,  and  the  use  of 
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such  methods  on  several  systems  of  chemical  interest.  Chapter  two  deals 
with  expanding  the  INDO  method  so  that  the  prediction  of  /-orbitals 
being  localized  orbitals  may  be  verified  and  so  that  the  unique  bonding 
and  spectroscopy  of  these  compounds  may  be  examined.  By  adapting  the 
INDO  method,  we  may  now  expand  our  studies  to  include  the  chemistry  of 
the  lanthanides  and  actinides. 

The  chemistry  of  the  lanthanides  and  actinides  is  different  from  the 
chemistry  of  the  corresponding  d-orbital  chemistry.  The  compact 
(localized)  nature  of  the  /-orbitals,  causes  the  /-/  spectral 
transitions  to  be  characterized  by  very  sharp  transitions  and  the 
positions  of  the  transitions  are  almost  unaffected  by  the  ligands 
attached  to  the  metal.   The  f-orbitals  are  potentially  involved  in 
expanding  the  valence  of  lanthanide  containing  compounds;  some 
lanthanide  molecules  have  a  coordination  number  of  nine  and  several 
twelve  coordinate  lanthanide  compounds  are  known.  Are  /-orbitals 
required  for  greater  valency,  or  is  the  greater  valency  merely  a 
consequence  of  the  larger  ionic  radius  of  most  lanthanides?  The  study 
by  quantum  chemical  methods  has  been  slowed  by  the  size  of  the 
lanthanide  containing  molecules,  but  the  INDO  method  lends  itself  to  the 
study  of  large  molecules  and  therefore  the  choice  was  made  to  expand  the 
INDO  model  to  include  /-orbitals. 


CHAPTER  ONE 
LOCALIZED  ORBITALS 

Background 

The  observable  properties  of  any  wavefunction  composed  of  a  single 
Slater  determinant  are  invariant  to  a  unitary  transformation  of  the 
orbitals  occupied  in  the  wavefunction.^  Because  of  this  invariance,  the 
observable  properties  of  a  closed-shell  self-consistent  field  (SCF) 
wavefunction  may  be  described  using  canonical  orbitals,  or  any  set  of 
orbitals  related  to  the  canonical  orbitals  by  a  unitary  transformation. 

Canonical  orbitals  are  quite  useful  in  post-Hartee-Fock  calculations 
for  several  reasons.   Canonical  molecular  orbitals  (CMOs)  are  obtained 
directly  by  matrix  diagonalization  from  the  SCF  procedure  itself.   The 
canonical  orbitals  form  irreducible  representations  of  the  molecular 
point  group.   Since  the  symmetry  is  maintained,  all  subsequent 
calculations  may  be  simplified  by  the  use  of  symmetry.   Spectroscopic 
selection  rules  are  determined  using  the  canonical  orbitals.   Koopman's 
theorem,  which  relates  orbital  energies  to  molecular  ionization 
potentials  and  electron  affinities,  is  based  entirely  on  the  use  of 
canonical  orbitals. 

Localized  orbitals  (LMOs)  allow  for  the  wavefunction  to  be 
interpreted  in  terms  of  bond  orbitals,  lone-pair  orbitals  and  inner- 
shell  orbitals,  consistent  with  the  Lewis  structures  learned  in  freshmen 
chemistry.   Unlike  CMOs,  LMOs  may  be  transferred  into  other 
wavefunctions  as  an  initial  gu.  ss,  thereby  reducing  the  effort  needed  to 
produce  wavefunctions  for  large  molecules.   The  most  important  use  of 
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localized  orbitals  is  their  ability  to  simplify  configuration 
interaction  (CI)  calculations.  LMOs  maximize  intra-orbital  electronic 
interactions  and  minimize  inter-orbital  electronic  interactions.   This 
concentrates  correlation  energy  into  several  large  portions  instead  of 
many  small  portions  as  given  by  the  CMOs .   A  major  disadvantage  of  the 
use  of  localized  orbitals  is  the  loss  of  molecular  point  group  symmetry. 
Localized  orbitals  do  not  transform  as  an  irreducible  representation  of 
the  molecular  point  group.   The  total  vavefunction,  of  course,  does. 

Localization  methods  may  be  divided  into  several  categories.  The 
first  category  of  localization  is  based  on  an  implicit  definition  of 
what  a  localized  orbital  should  be.   An  underlying  physical  basis  for 
localized  orbitals  is  exploited  in  the  second  category  of  localization. 
Localized  orbitals  may  also  be  produced  in  accord  with  the  users  own 
definition  of  localization. 

The  implicit  definition  on  which  the  localized  orbitals  are  produced 
differs  from  method  to  method  but  all  of  these  methods  proceed  in  a 
similar  fashion.   A  function  of  the  form 

n 

G  =  Z  <ii|g  |ii>,  (1_1) 

i  =  l     ii 

is  maximized  or  minimized,  where  the  definition  of  <ii|g..|ii>  depends 
on  the  localization  criterion.   One  choice  for  the  value  of 
^^■^ '^ijkl'^"'"''  ^^  '"^^  two-electron  repulsion  integrals;""  for  this 
choice  the  sum  G  is  maximized.   This  method  is  referred  to  as  the 
Edmiston  Ruedenberg  (ER)  method.   Perhaps  the  most  popular  choice  of  a 
localization  method  is  the  Boys  method,  in  which  the  g  operator  is  the 
orbital  self-extension  operator, 

^ii  =  ""U      •  (1-2) 
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This  form  of  the  g  operator  may  be  recast  in  terms  of  the  product  of  two 
molecular  orbital  dipole  operators.   One  can  relate  this  form  of 
localization  to  maximizing  the  distance  between  the  orbital  centroids. 

Once  a  localization  criterion  has  been  established,  the  next  step  is 
to  construct  a  transformation  matrix  to  do  the  localization.   Since  the 
exact  nature  of  the  transformation  is  unknown,  an  iterative  procedure  is 
used  to  construct  the  localized  orbitals.   This  iterative  procedure 
moves  from  a  less-localized  set  of  orbitals  to  a  more-localized  set. 
Once  a  convergence  criterion  is  met  i.e.,  the  orbitals  do  not  change 
within  a  given  tolerance,  the  iterative  procedure  is  stopped. 

Although  localizations  using  either  of  the  above  two  methods  are 
relatively  standard,  some  problems  may  be  encountered.   As  with  any 
iterative  procedure,  convergence  difficulties  may  be  encountered.   In 
the  case  of  the  ER  method  all  two-electron  repulsion  integrals  must  be 

transformed  on  each  iteration,  a  very  time  consuming  step  proportional 

5    .  -> 

to  N  .   Since  the  SCF  procedure  itself  proceeds  as  N   (semi-empirical) 

4 
or  N  (ab  initio)  and  the  systems  studied  here  are  large,  we  will  not 

consider  the  ER  method  of  localization  any  further.   The  same  integral 
transformation  problem  is  encountered  for  the  Boys  method,  but  since  the 
integrals  involved  are  dipole  (one-electron)  integrals  the  problem  is 
much  simpler.   Since  the  localization  criteria  are  so  different  there  is 
no  reason  to  expect  different  methods  to  yield  orbitals  that  are 
similar,  but  in  general  the  LMOs  are  quite  similar  for  the  Boys  and  ER 
methods.   These  orbital  similarities  lead  to  the  :;econd  catecrotrv  of 
localization. 

We  claim  that  the  undsrlyinc  physic.^J  has;.-  of  l^raLi^ation  is  the 
rermi  hoJp.      rhe  FiM.mi  hiil'^  pro-'idfv  a  di  vr- r  ( ii,.ri-i  tera  five)  method 
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for  transforming  canonical  orbitals  to  localized  orbitals.  The  integral 
transformations  that  limit  the  usefulness  of  the  Boys  and  ER  methods  are 
also  eliminated  when  using  this  method.  The  disadvantage  of  this  method 
is  the  fact  that  a  series  of  probe  points  must  be  generated  for  the 

molecule.   These  points  may  be  generated  using  chemical  intuition  or  by 

12 
a  search  of  the  Fermi  hole  mobility  function.    The  Fermi  hole  method 

of  localization  may  also  fall  into  the  final  category  since  it  can  be 

made  to  pick  out  a  particular  localized  orbital  set. 

The  final  category  of  localization  method  allows  one  to  produce 

orbitals  in  accordance  with  one's  needs.   As  mentioned  above,  the  Fermi 

method  may  be  classified  in  this  category,  but  another  method  was 

developed  especially  for  this  purpose,  one  that  we  have  called  the 

double  projector  (DP)  method.   This  method  has  been  used  in  conjunction 

with  the  other  methods  above  to  help  predict  the  lowest  energy  state  of 

lanthanide  containing  compounds  where  f-orbital  degeneracies  are  a 

problem.   The  DP  method  allows  one  to  separate  the  /-orbitals  from  the 

other  metal  orbitals  and  use  a  small  CI  to  determine  the  ground  state  of 

the  molecule. 

Double  Projector 


The  double  projector  (DP)  method  of  localization  is  an  extremely 
useful  method  for  localizing  orbitals  when  the  form  of  the  localized 
orbitals  is  known  or  suspected  in  advance.   For  example,  if  one  would 
like  to  study  n-n  transitions  in  a  molecule,  a  full  localization  need 
not  be  done,  the  double  projector  may  be  used  to  isolate  (localize)  the 
n-type  orbitals.  A  subsequent  small  singles  CI  may  then  be  used  to  study 
only  the  n-n  transitions  and  thereby  elucidate  the  n-n  spectra. 
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Another  example  involves  the  localization  of  the  d-orbitals  in  a 
transition  metal  complex.   Because  of  accidental  degeneracies  between 
metal  d-orbitals  and  ligand  molecular  orbitals  (MOs),  the  atomic  d- 
orbitals  may  be  spread  out  in  many  canonical  orbitals.   A  large  CI  is 
then  required  to  restore  the  localized  nature  of  the  d  type  molecular 
orbitals.   Such  a  large  CI  can  be  avoided  by  first  doing  a  DP 
localization. 

The  DP  method  can  also  be  used  to  remove  orbitals  from  the  orbital 
set  so  that  the  remaining  orbitals  may  be  localized  using  a  standard 
localization  technique.  For  example,  a  common  problem  with  a  Boys 
localization  is  the  mixing  of  a   and  n  orbitals  to  obtain  x  orbitals, 
this  is  not  desirable  since  the  a   and  ix  spectra  will  now  be  mixed  and 
more  difficult  to  interpret.   The  n   orbitals  may  be  removed  using  the  DP 
method,  the  remaining  orbitals  localized,  and  the  n-ii  spectra 
calculated  using  a  small  singles  CI.   The  double  projector  is  a 
complementary  method  of  localization  and  is  normally  used  in  conjunction 
with  other  traditional  methods  of  localization;  therefore,  no  examples 
of  its  use  will  be  given  here. 

An  outline  of  the  double  projector  method  is  given  in  this  section. 
Consider  a  set  of  m  occupied  spin  orbitals  {$.)^,  and  a  set  of  r 
localized  "pattern"  orbitals  {Y.}^,  where  r  is  less  than  or  equal  to  m. 
These  "pattern"  orbitals  are  projected  out  of  the  set  [^■],    by 

m 

|Y'>  =  Z      I^.Xf.  It  >  (1-3) 

'  a    .  1    1    1  '  a  ^    ' 

1  =  1 

J- 

for  a  =  1  to  r.   These  {|y>}^,  are  then  symmetrically  orthogonalized 


r  T  =  A 


e  =  T'  A" 


•1/2 


and  are  projected  out  of  the  original  set  {$.}'}', 


♦'  =  ♦  (  1  -  Z  19  ><0  I)   . 
a   a   a' 


The  matrix  A'  is  formed  and  diagonalized 
U  i^i'\J   =  U'^A  U  =  X 


(1-4) 
(1-5) 

(1-6) 


(1-7) 


The  X  matrix  will  have  r  near  zero  eigenvalues  corresponding  to  the 
^®a^l  ^^^^  ^^^®  ^^^"  projected  out.   These  eigenvalues  and  the 
corresponding  columns  of  U  are  removed.   The  new  set  of  orthonormalized 


,  m-r 


orbitals  {Y  ]'    is  formed  from 
a  1 


♦'UX' 


•1/2 


(1-8) 


This  set  is  an  orthogonal  complement  to  the  set  10  >,  but  has  no 

a 

particular  physical  significance.  To  obtain  a  set  of  orbitals  most  like 
the  canonical  set,  we  form  F,  the  Fock  matrix,  over  the  Y  subset  and 
diagonalize  F, 


W'^Y'^FYW  =  e"^ 


Y'  =  YW 


(1-9) 
(1-10) 


Y'  are  linear  combinations  of  Y  that  we  can  energy  order  according  to 
e  .   These  Y'  orbitals  are  the  most  like  the  original  canonical  orbitals 
with  the  "pattern"  orbitals  removed. 


t".. 
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Fermi  Localization 

Background 

This  section  presents  a  method  for  transforming  a  set  of  canonical 
SCF  orbitals  into  a  set  of  localized  orbitals  based  on  the  properties  of 
the  Fermi  hole  ~   and  the  Fermi  orbital.  ■'■■^' ■'"'^  Unlike  localization 
methods  based  on  iterative  optimization  of  some  criterion  of 
localization,    '   '    the  method  presented  here  provides  a  direct  (non- 
iterative)  calculation  of  the  localized  orbital  transformation  matrix. 
Consequently,  this  method  avoids  the  convergence  problems  which  are 
possible  with  iterative  transformations. 

Unlike  the  extrinsic  methods  for  transforming  canonical  SCF  orbitals 

17-19 
into  localized  orbitals,      the  method  presented  here  does  not  depend 

on  the  introduction  of  a  definition  of  a  set  of  "atomic  orbitals".   The 

method  presented  here  may  also  be  distinguished  from  applications  of 

localized  orbitals  such  as  the  PCLIO  method      in  that  the  latter 

method  does  not  involve  SCF  orbitals,  and  it  is  nor  concerned  with  the 

transformation  of  canonical  SCF  orbitals  into  localized  orbitals. 

Properties  of  the  Fermi  Hole 

The  Fermi  hole  is  defined  as 

A(rj;r2)  =  p(r^)  -  2p^{r^,r^)/ pit^) ,  (1-11) 

where  p(r^)  is  the  diagonal  portion  of  the  first  order  reduced  density 
matrix  and  p(r^;r2)  is  the  corresponding  part  of  the  second  order 
reduced  density  matrix.    For  special  case  of  a  closed  shell  SCF 
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vavefunction,  the  natural  representation  of  the  Fermi  hole  is  the 

IT  1  / 

absolute  square  of  the  Fermi  orbital   ' 

2 


The  Fermi  orbital  is  given  by 

f(rT;r„)  =  [2/p(r„)]^^2  Z  g.  (r^  )g.  (rO  , 


(1-12) 


•l'"2 


i'  I'^i^  2' 


(1-13) 


where  the  orbitals  g.(r)  are  either  the  canonical  SCF  molecular  orbitals 
or  any  set  related  to  the  canonical  SCF  molecular  orbitals  by  a  unitary 
transformation.   The  Fermi  orbital  f(r^;r„)  is  interpreted  as  a  function 
of  r^  which  is  parametrically  dependent  upon  the  position  of  a  probe 

electron  l(icated  at  x^. 

12  13  27  28 
Previous  work   '   '   '   has  demonstrated  that  the  Fermi  hole  does 

not  follow  the  probe  electron  in  a  uniform  manner.   Instead,  molecules 

are  found  to  possess  regions  where  the  Fermi  hole  is  insensitive  to  the 

position  of  the  probe  electron.   As  the  probe  electron  passes  through 

one  of  these  regions,  the  Fermi  hole  remains  nearly  stationary  with 

respect  to  the  nuclei.   These  regions  are  separated  by  regions  where  the 

Fermi  hole  is  very  sensitive  to  the  position  of  the  probe  electron.   As 

the  probe  electron  passes  through  one  of  these  regions,  the  Fermi  hole 

changes  rapidly  from  one  stable  form  to  another. 

The  sensitivity  of  the  Fermi  hole  to  the  position  of  the  probe 

electron  is  measured  by  the  Fermi  hole  mobility  function,   '' 


where 


F(r)  =  F^(r)  +  Fy(r)  +  F^(r) 


2 


(1-U) 


^(•^)  =   2 
P 


> 

^i 

3vJ 

-^J 

3v 

i>j 

(1-15) 


for  V  =  X,  y  or  z.   This  may  be  compared  to 
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Fq(p)  =  (3ti/4)(p/2)^^^  (1-16) 

which  provides  an  estimate  of  the  Fermi  hole  in  a  uniform  density 
electron  gas. 

The  Fermi  hole  mobility  function  F(r)  for  the  formaldehyde  molecule 
is  shown  in  Fig.  1-1.   The  difference  F(r)-F^(p)  is  shown  in  Fig.  1-2. 
Regions  where  F(r)  >  F(p)  ,  that  is,  the  Fermi  hole  is  less  sensitive  to 
the  position  of  the  probe  electron  than  it  would  be  in  an  electron  gas 

of  the  same  density,  may  be  compared  to  the  loges  proposed  by 

29-33 
Daudel.       Regions  where  F(r)  =  F(p)  resemble  boundaries  between 

loges. 

When  the  probe  electron  is  located  in  a  region  where  F(r)  <  F(p), 

the  Fermi  orbital  is  found  to  resemble  a  localized  orbital  determined  by 

conventional  methods.  "  '   '    This  similarity  is  demonstrated  by  Figs. 

1-3  and  1-4  which  compare  a  Fermi  hole  for  the  formaldehyde  molecule 

with  a  localized  orbital  determined  by  the  orbital  centroid  criterion  of 

-.        T   ,.   4-6,15 
localization.    ' 

Localized  Orbitals  Based  on  the  Fermi  hole 


Equation  1-13  provides  a  direct  relationship  between  a  set  of 
canonical  SCF  orbitals  g.(r)  and  a  localized  orbital  f.(r)  =  f(r,r.) 
where  r.  is  a  point  in  a  region  where  F(r.)  <  F^(p(r.))-   In  order  to 
transform  a  set  of  N  canonical  SCF  orbitals  into  a  set  of  N  localized 
orbitals,  it  is  necessary  to  select  N  points  r.  j  =  1  to  N,  each  of 
which  is  located  in  a  region  where  F(r.)  <  F„(p(r.))-   Ideally,  each  of 
these  points  should  correspond  to  a  minimum  of  F(r)  or  F(r)-Fp.(p).   This 
condition,  however,  is  not  critical,  because  the  Fermi  hole  is 
relatively  insensitive  to  the  position  of  the  probe  electron  when  the 
probe  electron  is  located  in  one  of  these  regions. 


t'T- 
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Figure  1-1:  The  fermi  hole  mobility  function  F(r)  for  the  H„CO  based  on 
the  geometry  and  double  zeta  basis  set  of  ref.  41.   The 
locations  of  the  nuclei  are  indicated  by  (+)  signs.   The 
contours  represent  mobility  function  values  of  0.1,  0.25, 
0.5,  1.0,  2.0  and  5.0  atomic  units.   The  contours  increase 
from  0.1  near  the  corners,  to  over  5.0  in  regions  enclosing 
the  carbon  and  oxygen  nuclei.   Each  nucleus  is  located  at  a 
local  minimum  of  the  mobility  function. 


*A 
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Figure  1-2:  The  difference  between  the  Fermi  hole  mobility  function  F(r) 
and  the  electron  gas  approximation  for  the  H„CO  molecule. 
The  contours  represent  values  of  0.0,  -0.1,  -0.25,  -0.5, 
-1.0,  -2.0  and  -5.0,  in  addition  to  those  indicated  in 
figure  1-1.   the  contours  representing  negative  values  and 
zero  are  indicated  by  broken  lines.   Each  nucleus  is  located 
at  a  local  minimum. 


u 


Figure  1-3:  The  fermi  hole  for  the  formaldehyde  molecule  determined  by  a 
probe  electron  located  at  one  of  the  protons.   The  contours 
indicate  electron  density  of  0.005,  0.01,  0.02,  0.04,  0.08, 
0.16,  0.32,  0.6A,  1.28  and  2.56  electrons  per  cubic  bohr. 
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Figure  1-4:  The  localized  orbital  for  the  C-H  bond  of  a  formaldehyde 
molecule  determined  by  the  orbital  centriod  criterion  for 
localization.   The  electronic  density  contours  are  the  same 
as  in  figure  1-3. 
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A  set  of  N  Fermi  orbitals  determined  by  Eq.  1-3  is  not  generally 
orthonormal.   Each  member  of  this  set,  however,  is  usually  very  similar 
to  one  member  of  an  orthonormal  set  of  conventional  localized  orbitals. 
Consequently,  the  overlap  between  a  pair  of  Fermi  orbitals  is  usually 
very  small,  and  a  set  of  N  Fermi  orbitals  may  easily  be  converted  into 

an  orthonormal  set  of  localized  orbitals  by  means  of  the  method  of 

33 
symmetric  orthogonalization.    The  resulting  unitary  transformation  is 

given  by 

U  =  (TT^)"^^^T,  (1-17) 

where 

T^j  =  gj(r.)/(p(r.)/2)^^^  (1-18) 

In  the  following  three  sections,  the  transformation  of  canonical  SCF 
orbitals  based  on  Eqs.  1-17  and  1-18  is  demonstrated  for  each  of  three 
molecules.  The  first  example,  a  cyclic  conjugated  enone,  represents  a 
simple  case  where  conventional  methods  are  not  expected  to  have  any 
special  difficulties.   The  second  example,  methyl  acetylene,  is  a 

molecule  for  which  conventional  methods  have  serious  convergence 

34 
problems.    The  third  example,  boron  trifluoride,  is  a  pathological 

case  for  the  orbital  centroid  criterion,  with  a  number  of  local  maxima 

,  and  saddle  points  in  the  potential  surface  according  to  the  Boys 

criterion  of  localization. 

In  each  case,  the  first  step  in  the  application  of  this  method  is 

the  selection  of  the  set  N  points.   This  set  always  includes  the 

locations  of  all  of  the  nuclei  in  the  molecule.   For  atoms  other  than 

hydrogen,  the  resulting  Fermi  orbitals  are  similar  to  innershell 

localized  orbital?.   When  the  probe  electron  is  located  on  a  hydrogen 

l]  atom,  the  Fermi  orbital  is  similar  to  an  R-H  bond  orbital. 
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Additional  points  for  the  probe  electron  may  usually  be  determined 
based  on  the  molecular  geometry.   The  midpoint  between  two  bonded  atoms 
(other  than  hydrogen)  tends  to  yield  a  Fermi  orbital  resembling  a  single 
bond.   Multiple  bonds  may  be  represented  with  two  or  three  points 
located  roughly  one  to  two  bohr  from  a  point  midway  between  the  multiply 
bonded  atoms,  along  lines  perpendicular  to  a  line  joining  the  nuclei. 
Likewise,  lone  pair  orbitals  may  be  determined  by  points  located  roughly 
one  bohr  from  the  nucleus  of  an  atom  which  is  expected  to  possess  lone 
pair  orbitals. 

Application  to  the  Furanone  Molecule 

The  furanone  molecule,  C.H.O^,  and  its  derivatives  are  useful 

35-37 
reagents  in  2+2  photochemical  cycloadditions.       The  canonical  SCF 

molecular  orbital  for  the  furanone  molecule  were  calculated  with  an  STO- 

38 

3G   basis  set  and  the  geometry  specified  in  Table  1-1.   The  molecular 

geometry  was  restricted  to  C  symmetry,  with  a  planar  five  membered 
ring.   Fermi  hole  localized  orbitals  were  calculated  based  on  the  set  of 
points  indicated  in  Table  1-1.   These  points  include  the  positions  of 
the  ten  nuclei,  as  well  as  twelve  additional  points  determined  by  the 
method  outlined  above. 

The  centroids  of  the  localized  orbitals  determined  by  the  points  in 
Table  1-1  are  shown  in  Table  1-2.   The  C=C  and  C=0  double  bonds  are  each 
represented  by  a  pair  of  equivalent  bent  (banana)  bonds  similar  to  those 
determined  by  other  methods  for  transforming  canonical  SCF  orbitals  into 
localized  orbitals. 

As  shown  in  Table  1-2,  the  centroids  of  the  orbitals  determined  by 
the  Fermi  hole  are  very  close  to  those  of  the  localized  orbitals  Page 
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Table  1-1  :  Molecular  geometry  and  probe  electron  points  for  the 
furanone  (C.H.O^).   The  first  ten  points  indicate  the 
molecular  geometry  used  in  these  calculations.   The  twelve 
additional  probe  electron  positions  were  determined  as 
described  in  the  text.   All  coordinates  are  given  in  bohr. 


Position 

X 

Y 

Z 

Atom  C^ 

0.0 

0.0 

0.0 

Atom  C2 

0.0 

2.589 

0.0 

Atom  C^ 

2.671 

3.355 

0.0 

Atom  C, 
4 

4.363 

0.940 

0.0 

Atom  0, 

2.534 

-1.071 

0.0 

Atom  Op 

3.427 

5.548 

0.0 

Atom  H^ 

-1.802 

-1.040 

0.0 

Atom  H„ 

-1.802 

3.629 

0.0 

Atom  H-j 

5.560 

0.858 

1.698 

Atom  H, 

4 

5.560 

0.858 

-1.698 

C.-C2  bond  1 

0.0 

1.295 

2.000 

C,-C2  bond  2 

0.0 

1.295 

-2.000 

C2-C2  bond 

1.336 

2.972 

0.0 

C^-C^  bond 

3.517 

2.148 

0.0 

C,-0,  bond 
4   1 

3.449 

-0.066 

0.0 

Cj-0,  bond 

1.267 

-0.536 

0.0 

CyO^   bond  1 

3.049 

4.452 

2.000 

0^-02  bond  2 

3.049 

4.452 

-2.000 

0,  lone  pair  1 

2.654 

-1.821 

0.660 

0.  lone  pair  2 

2.654 

-1.821 

-0.660 

Oj   lone  pair  1 

2.707 

6.298 

0.0 

Op  lone  pair  2 

4.267 

5.698 

0.0 
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Table  1-2  :  Orbital  centroids  for  localized  orbitals  determined  by  the 

Fermi  hole  method  and  by  the  orbital  centroid  method  for  the 
furanone  molecule  (C,H,0„).  All  coordinates  are  given  in 
bohr. 


Orbital 

Fermi  hole  method 

Centroid  criterion 

X 

Y 

Z 

X 

Y 

Z 

C^  K  shell 

0.0 

0.001 

0.0 

0.0 

0.0 

0.0 

C^   K  shell 

0.0 

2.588 

0.0 

0.0 

2.588 

0.0 

C^  K  shell 

2.670 

3.355 

U.O 

2.671 

3.355 

0.0 

C^  K  shell 

4.362 

0.940 

0.0 

4.326 

0.939 

0.0 

0^  K  shell 

2.534 

-1.070 

0.0 

2.533 

-1.070 

0.0 

O2  K  shell 

3.427 

5 .  547 

0.0 

3.426 

5.547 

0.0 

C^-H^  bond 

-1.138 

-0.742 

0.0 

-1.171 

-0.761 

0.0 

C2-H2  bond 

-1.148 

3.318 

0.0 

-1.181 

3.331 

0.0 

C,-Ht  bond 
4  3 

5.149 

0.895 

1.132 

5.155 

0.888 

1.153 

C,-H,  bond 
4   4 

5.149 

0.895 

-1.132 

5.155 

0.888 

-1.153 

C^-C^   bond  1 

0.008 

1.401 

0.635 

0.030 

1.410 

0.599 

C.-C^  bond  2 

0.008 

1.401 

-0.635 

0.030 

1.410 

-0.599 

C2-C^  bond 

1.320 

3.054 

0.0 

1.284 

3.046 

0.0 

C^-C^  bond 

3.542 

2.190 

0.0 

3.565 

2.149 

0.0 

C.-O,  bond 
4  1 

3.268 

-0.243 

0.0 

3.243 

-0.225 

0.0 

C^-0^  bond 

1.472 

-0.627 

0.0 

1.495 

-0.613 

0.0 

0^-02  bond  1 

3.093 

4.584 

0.548 

3.114 

4.612 

0.511 

C2-O2  bond  2 

3.093 

4.584 

-0.548 

3.114 

4.612 

-0.511 

0.  lone  pair  1 

2.554 

-1.283 

0.440 

2.594 

-l.?il 

0.472 

0^  lone  pair  2 

2.554 

-1.283 

-0.440 

2.594 

-1.311 

-0.472- 

O2  lone  pair  1 

3.015 

5.886 

0.0 

3.011 

5.892 

0.0 

O2  lone  pair  2 

3.970 

5.589 

0.0 

3.970 

5.600 

0.0 
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determined  by  the  orbital  centroid  criterion.    '    Likewise,  the 
localized  orbitals  determined  by  the  Fermi  hole  were  found  to  be  very 
close  to  those  determined  by  the  orbital  centroid  criterion.   Each  of 
the  localized  orbitals  determined  by  the  Fermi  hole  method  was  found  to 
have  an  overlap  of  0.994  to  0.999  with  one  of  the  localized  orbitals 
determined  by  the  orbital  centroid  criterion.  The  remaining  (off- 
diagonal)  overlap  integrals  between  these  two  sets  of  localized  orbitals 
were  found  to  have  a  root  mean  square  (RMS)  value  of  0.011734. 

The  transformation  of  a  set  of  canonical  SCF  orbitals  to  an 
orthonormal  set  of  localized  orbitals  determined  by  the  Fermi  hole 
required  10  minutes  on  a  PDP-11/44  computer.   The  orbital  centroid 
(Boys)  method  required  140  minutes  starting  from  the  canonical  SCF 
molecular  orbitals  or  80  minutes,  using  the  Fermi  localized  orbitals  as 
an  initial  guess,  to  reach  T^^^^^  of  less  than  10"^,  where  T^^j^  is  the  RMS 
value  of  the  oCf-diagonal  part  of  the  transformation  matrix  which 
converts  the  orbitals  obtained  on  one  iteration  to  those  of  the  next 
iteration.   The  orbital  centroid  criterion  calculations  reported  here 
are  based  on  a  partially  quadratic  procedure   which  requires  less  time 
than  conventional  localization  procedures  based  on  2X2  rotation. 

Application  to  Methylacetylene 


The  localized  orbitals  of  methylacetylene  are  of  interest  because  of 
the  convergence  difficulties  encountered  in  attempts  to  calculate  these 
orbitals  using  iterative  localization  methods.   These  difficulties  are 
caused  by  the  weak  dependence  of  the  criterion  of  localization  on  the 
orientation  of  the  three  equivalent  C-C  (banana)  bonds  relative  to  the 
i  '•        three  C-H  bonds  of  the  methyl  group.   In  calculations  based  on  the 
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orbital  centroid  criterion,  over  200  iterations  were  required  to 

determine  a  set  of  orbitals  which  satisfied  a  very  weak  criterion  of 

34 
convergence.    Most  of  these  difficulties  may  be  overcome  using  the 

quadratically  convergent  method  which  has  been  developed  recently. '''^  As 
shown  below,  however,  the  localized  orbitals  based  on  the  Fermi  hole 
yield  nearly  equivalent  results  and  require  much  less  effort  than  even 
the  quadratically  convergent  method. 

The  canonical  SCF  molecular  orbitals  for  methylacetylene  were 
determined  by  an  ST0-5G  basis  set^^^  and  an  experimental  geometry.  "^^ 
Transformation  of  the  11  occupied  SCF  orbitals  into  localized  orbitals 
based  on  the  Fermi  hole  method  required  the  selection  of  11  points. 
These  points  are  shown  in  Table  1-3.   The  positions  of  the  nuclei 
provided  seven  of  these  points.   One  point  was  located  at  the  midpoint 
of  the  C2-C^  single  bond.  The  remaining  three  points  were  located  two 
bohr  from  the  l>^  rotation  axis  at  a  point  midway  between  the  C-,  and  C„ 
nuclei.   These  last  three  points  were  eclipsed  with  respect  to  the 
methyl  protons. 

The  centroids  of  the  localized  orbitals  determined  by  this  method 
are  shown  in  Table  1-4.   As  expected,  the  triple  bond  is  represented 
with  three  equivalent  banana  bonds.  The  centroids  of  the  corresponding 
orbitals  determined  by  the  orbital  centroid  criterion  are  also  shown  in 
Table  1-4.   These  are  very  close  to  those  determined  by  the  Fermi  hole 
method.   The  RMS  value  of  the  off-diagonal  part  of  the  overlap  matrix 
between  the  localized  orbitals  determined  by  the  Fermi  hole  and  those 
determined  by  the  orbital  centroid  criterion  is  0.012874. 

The  Fermi  hole  method  required  1.62  minutes  to  transform  the 
,.>,       canonical  SCF  molecular  orbitals  into  an  orthonormal  set  of  localized 
:'M^      "lolecular  orbitals.   By  comparison,  the  (quadratically  convergent)  Page 
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Molecular  geometry  and  probe  electron  positions  for  the 
methylacetylene  molecule.   The  first  seven  points  indicate 
the  locations  of  the  nuclei.   All  coordinates  are  given  in 
bohr. 


Position 

X 

Y 

Z 

Atom  C. 

0.0 

-1.140 

0.0 

Atom  C2 

0.0 

1.140 

0.0 

Atom  C^ 

0.0 

3.897 

0.0 

Atom  H^ 

0.0 

-3.145 

0.0 

Atom  H- 

-1.961 

4.438 

0.0 

Atom  H- 

0.980 

4.438 

1.698 

Atom  H, 
4 

0.980 

4.438 

-1.698 

CJ-C2  bond  1 

-2.0 

0.0 

0.0 

C^-C2  bond  2 

1.0 

0.0 

1.732 

Cj-C2  bond  3 

1.0 

0.0 

-1.732 

C2-C2  bond 

0.0 

2.510 

0.0 

Table  1-4 
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Orbital  centroids  for  localized  orbitals  determined  by  the 
Fermi  hole  method  and  by  the  orbital  centroi-l  method  for 
the  methylacetylene  molecule.   All  coordinates  are  given 
in  bohr. 


Orbital 

Fermi  hole  method 

Centroid  criterion 

X 

Y 

Z 

X 

Y 

Z 

Cj  K  shell 

0.0 

-1.141 

0.0 

0.0 

-1.140 

0.0 

C2  K  shell 

0.0 

1.139 

0.0 

0.0 

1.134 

0.0 

C3  K  shell 

0.0 

3.896 

0.0 

0.0 

3.896 

0.0 

C^-Hj  bond 

0.0 

-2.486 

0.0 

0.0 

-2.502 

0.0 

Cy^2   ''°"^ 

-1.302 

4.270 

0.0 

-3.316 

4.272 

0.0 

C^-H^  bond 

0.651 

4.270 

1.127 

0.658 

4.272- 

1.139 

C^-H^  bond 

0.651 

4.270 

-1.127 

0.658 

4.272 

-1.139 

C^-C2  bond  1 

-0.713 

-0.003 

0.0 

-0.692 

-0.U13 

0.0 

Cj-C2  bond  2 

0.356 

-0.003 

0.617 

0.346 

-0.013 

0.599 

Cj-Cj  bond  3 

0.356 

-0.003 

-0.617 

0.346 

-0.013 

-0.599 

C»-C-)  bond 

0.0 

2.457 

0.0 

0.0 

2.498 

0.0 

I:- 

Si'': 
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orbital  centroid  method  required  26  minutes  starting  from  the  canonical 
SCF  molecular  orbitals  or  19  minutes  using  the  Fermi  localized  orbitals 
as  an  initial  guess  to  reach  T^^^   of  10"^  or  less. 

Application  to  Boron  Trifluoride 

As  a  further  example  of  the  application  of  the  Fermi  hole 
localization  method,  localized  orbitals  were  also  calculated  for  the 
boron  trifluoride  molecule.   The  molecule  provides  a  demonstration  of 
how  characteristics  of  little  or  no  physical  significance  can  cause 
serious  convergence  difficulties  for  iterative  localization  methods.   A 
localized  representation  of  boron  trifluoride  includes  four  inner-shell 
orbitals,  three  boron  fluorine  bond  orbitals,  and  nine  fluorine  lone 
pair  orbitals.   The  orbital  centroid  criterion  method  shows  a  small 
dependence  on  rotation  of  each  set  of  three  lone  pair  orbitals  about  the 
corresponding  B-F  axis.   Consequently,  the  hessian  matrix  for  he 
criterion  of  localization  as  a  function  of  a  unitary  transformation  of 
the  orbitals  has  three  very  small  eigenvalues. 

The  optimal  orientation  of  the  fluorine  lone  pairs  may  correspond  to 
one  of  several  possible  conformations.  One  of  these,  the  "pinwheel" 
conformation,  has  a  single  lone  pair  orbital  on  one  of  the  fluorine 
atoms  in  the  plane  of  the  molecule.   The  other  two  lone  pair  orbitals  on 
this  fluorine  atom  are  related  to  the  first  lone  pair  by  120  degree 
rotations  about  the  F-B  axis.   The  lone  pair  orbitals  on  the  other  two 
fluorines  are  obtained  by  120  degree  rotations  about  the  C,  axis.   The 

point  group  of  the  orbital  centroids  for  this  conformation  is  C 

3h 

The  "three-up"  conformation  is  generated  by  rotating  the  set  of  lone 
pair  orbitals  on  each  fluorine  atom  in  the  pinwheel  conformation  by  90 
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degrees  about  each  F-B  axis.   The  point  group  for  the  orbital  centroids 

of  this  conformation  is  C^^.  The  "up-up-down"  conformation  is  generated 

by  rotating  the  set  of  lone  pair  orbitals  on  one  of  the  fluorine  atoms 

in  the  "three-up"  conformation  by  180  degrees  about  the  F-B  axis.   This 

conformation  has  the  symmetry  of  the  C  point  group. 

The  canonical  SCF  orbitals  for  BF^  were  calculated  based  on  the 

double-zeta  basis  set  and  geometry  tabulated  by  Snyder  and  Basch.^° 

Localized  orbitals  determined  by  the  orbital  centroid  criterion  were 

obtained  for  the  three-up  conformation  and  the  up-up-down  conformation. 

The  centroids  for  these  orbitals  are  shown  in  Table  1-5.   The  first  five 

(most  positive)  eigenvalues  of  the  hessian  matrix  for  each  of  these 

conformations  are  shown  in  Table  1-6.   All  eigenvalues  of  the  hessian 

matrix  are  negative  for  both  of  these  conformations,  indicating  that 

both  conformations  are  maxima  for  the  sum  of  squares  of  the  orbital 

centroids.   The  pinwheel  conformation,  however,  was  never  found. 

Consequently,  it  was  not  possible  to  exclude  the  possibility  that  the 

pinwheel  conformation  was  the  global  maximum  and  the  three-up 

conformation  was  only  a  local  maximum. 

The  pinwheel  conformation  can  easily  be  constructed  using  the  Fermi 

hole  localization  method  by  selecting  an  appropriate  set  of  probe 

positions.   This  set  of  points  is  shown  in  Table  1-7.   The  centroids  of 

the  resulting  set  of  localized  orbitals  are  shown  in  Table  1-8.   When 

this  set  of  localized  orbitals  is  used  as  the  starting  point,  the 

orbital  centroid  method  quickly  converges  to  a  stationary  point  with  C^^ 

3h 

symmetry.   The  centroids  of  the  resulting  set  of  orbitals  are  shown  in 
Table  1-8.   As  shown  in  Table  1-6,  three  of  the  eigenvalues  of  the 
hessian  matrix  were  positive  at  this  point,  demonstrating  that  the  Page 


^:, 
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Table  1-5 
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Orbital  centroids  for  localized  orbitals  determined  by  the 
orbital  centroid  method  for  the  boron  trifluoride 
molecule.   The  four  innershell  orbitals  have  been  excluded 
from  these  calculations.   All  coordinates  are  given  in 
bohr. 


Orbital 

up 

-up-up  (C 

3v) 

up- 

up-down  (C  ) 

0 

X 

Y 

Z 

X 

Y 

Z 

B-F^  bond 

1.725 

0.0 

-0.072 

1.725 

0.0 

-0.070 

B-F2  bond 

-0.862 

1.494 

-0.072 

-0.862 

1.494 

-0.070 

B-F^  bond 

-0.862 

-1.494 

-0.072 

-0.862 

-1.494 

0.069 

F^  lone  pair  1 

2.445 

0.0 

0.487 

2.446 

0.0 

0.487 

lone  pair  2 

2.596 

0.433 

-0.207 

2.596 

0.433 

-0.208 

lone  pair  3 

2.596 

-0.433 

-0.207 

2.596 

-0.433 

-0.207 

F2  lone  pair  1 

-1.222 

2.117 

0.487 

-1.222 

2.118 

0.487 

lone  pair  2 

-0.922 

2.465 

-0.207 

-0.922 

2.465 

-0.208 

lone  pair  3 

-1.674 

2.032 

-0.207 

-1.674 

2.031 

0.207 

F„  lone  pair  1 

-1.222 

-2.117 

0.487 

-1.223 

-2.119 

-0.486 

lone  pair  2 

-0.922 

-2.465 

-0.207 

-0.922 

-2.464 

0.208 

lone  pair  3 

-1.674 

-2.032 

-0.207 



-1.673 

-2.031 

0.208 
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Table  1-6  :  Values  of  the  orbital  centroid  criterion  and  the  second 

derivatives  of  the  orbital  centroid  criterion  for  various 
conformations  of  localized  orbitals  for  the  boron 
trifluoride  molecule.  The  row  labelled  sum  indicated  the 
sura  of  the  squares  of  the  orbital  centriods  for  each  of 
the  conformations.   The  following  rows  show  the  five 
highest  (most  positive)  eigenvalues  X.  of  the 
corresponding  hessian  matrix.  The  first  and  second 
columns  correspond  to  the  localized  orbitals  described  in 
described  in  Table  1-5.   The  third  and  fourth  columns 
correspond  to  localized  orbitals  described  in  Table  1-8. 
The  gradient  vectors  are  zero  for  the  first  three  columns. 


Configurations 

up-up-up 

up-up  down 

pinwheel 

Fermi  hole 

Sum 

69.445360 

69.444985 

69.438724 

69.341094 

\ 

-0.019521 

-0.018174 

+0.017234 

+0.031444 

^2 

-0.019850 

-0.018914 

+0.016238 

+0.029932 

h 

-0.019850 

-0.019095 

+0.016236 

+0.029928 

\ 

-0.155480 

-0.153827 

-0.193521 

-0.195502 

h 

-0.155484 

-0.155181 

-0.1'^4366 

-0.195524 

m 
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Table  1-7 


:  I'robe  electron  positions  for  the  boron  trifluoride 
molecule.   The  first  three  points  are  located  at  the 
midpoint  of  the  B-F  bonds.   The  remaining  points  have  been 
chosen  in  the  pinvheel  conformation  (symmetry  C^,  ).   All 
coordinates  are  given  in  bohr. 


Position 

X 

Y 

Z 

B-F^  bond 

1.223 

0.0 

0.0 

B-F,  bond 

-0.611 

1.059 

0.0 

B-F^  bond 

-0.611 

-1.059 

0.0 

F.  lone  pair  1 

2.781 

0.943 

0.0 

lone  pair  2 

2.718 

-0.472 

0.817 

lone  pair  3 

2.781 

-0.472 

-0.817 

F-  lone  pair  1 

-2.207 

1.937 

0.0 

lone  pair  2 

-0.981 

2.644 

0.817 

lone  pair  3 

-0.981 

2.644 

-0.817 

F-  lone  pair  1 

-0.573 

-2.880 

0.0 

lone  pair  2 

-1.799 

-2.172 

0.817 

lone  pair  3 

-1.799 

-2.172 

-0.817 

^; 


:-  I   -•  u. 


Table  1-8 
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Orbital  centroids  for  localized  orbitals  determined  by  the 
Fermi  hole  method  and  by  the  orbital  cr-ntroid  method  for 
the  boron  trJfluoiidt:  molecule.   The  tuur  innershell 
orbitals  have  been  excluded  from  these  calculations.   All 
coordinates  are  given  in  bohr. 


Orbital 

Fermi  hole 

Centroid  criterion 

X 

Y 

Z 

X 

Y 

Z 

B-F^  bond 

1.713 

0.0 

0.0 

1.720 

O.028 

0.0 

B-F^  bond 

-0.857 

1.484 

0.0 

-0.884 

1.475 

0.0 

B-F^  bond 

-0.837 

1.4n4 

0.0 

-0.835 

1.504 

0.0 

F-  lone  pair  1 

2.567 

0.488 

0.0 

2.603 

0.484 

0.0 

lone  pair  2 

2.538 

-0.244 

0.414 

2.520 

-0.256 

0.410 

lone  pair  3 

2.538 

-0.244 

-0.414 

2.520 

-0.2''6 

-0.410 

^2   T-one  pair  1 

-1.706 

1.976 

0.0 

-1.721 

2.012 

0.0 

lone  pair  2 

-1.057 

2.320 

0.414 

-1.037 

2.310 

0.410 

lone  pair  3 

-1.057 

2.32U 

-0.414 

-1.037 

2.310 

-0.410 

F~  lone  pair  1 

-0.861 

-2.467 

0.0 

-0.882 

-2.497 

0.0 

lone  pair  2 

-1.480 

-2.076 

0.414 

-1.482 

-2.054 

0.410 

lone  pair  3 

-1.480 

- 

-2.076 

-0.414 

-1.482 

-2.054 

-0.410 
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pinwheel  conformation  is  a  saddle  point  with  respect  to  the  orbital 
centroid  criterion.   These  calculations  also  indicate  that  the  three-up 
conformation  is  probably  the  global  maximum  for  the  orbital  centroid 
criterion. 

We  do  not  intend  to  attribute  any  special  physical  significance  to 
any  of  the  lone  pair  configurations  for  BF^.   These  calculations 
demonstrate  some  of  the  problems,  such  as  local  maxima  and  saddle 
points,  which  may  occur  for  conventional  iterative  localization  methods. 
These  calculations  demonstrate  how  the  Fermi  hole  method  may  be  used  by 
itself  to  transform  the  canonical  SCF  oibitals  into  localized  orbitals 
without  any  of  these  difficulties.  In  addition,  these  calculations 
demonstrate  how  the  Fermi  hole  method  may  be  used  in  conjunction  with 
the  orbital  centroid  method  to  establish  a  characteristic  of  the  orbital 
centroid  criterion  which  would  have  been  very  difficult  to  establish 
using  the  orbital  centroid  method  alone. 

Conclusions 

The  numerical  results  presented  here  demonstrate  how  the  properties 
of  the  Fermi  hole  may  be  used  to  transform  canonical  SCF  molecular 
orbitals  into  a  set  of  localized  SCF  molecular  orbitals.  Except  for  the 
symmetric  orthogonalization,  this  method  requires  no  integrals  and  no 
iterative  transformations.   The  localized  orbitals  obtained  from  this 
method  are  very  similar  to  the  localized  orbitals  determined  by  the 
orbital  centroid  criterion.   The  orbitals  determined  by  the  Fermi  hole 
may  be  used  directly  in  subsequent  calculations  requiring  localized 
orbitals.   Alternatively,  the  orbitals  determined  by  this  method  may  be 
used  as  a  starting  point  for  iterative  localization  procedures. ^"^' ^^' ^^ 
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The  necessity  of  providing  the  set  of  probe  electron  positions  may 
appear  to  introduce  a  subjective  element  into  the  localized  orbitals 
determined  by  the  Fermi  hole  method.   Most  of  the  subjective  character 
to  this  choice,  hovever,  is  eliminated  by  the  fact  that  Fermi  hole  is 
relatively  insensitive  to  the  location  of  the  probe  electron  whenever 
the  probe  electron  is  located  in  a  region  associated  with  a  strongly 
localized  orbital.   This  is  reflected  by  the  fact  that  the  centroids  of 
the  localized  orbitals  determined  by  the  Fermi  hole  method,  as  shown  in 
Tables  1-2,  1-4  and  1-8,  are  much  closer  to  the  centroids  of  the 
corresponding  localized  orbitals  determined  by  the  orbital  centroid 
criterion  than  they  are  to  the  probe  electron  positions  used  to 


calculate  them. 


If  the  electrons  are  not  strongly  localized  in  certain  portions  of  a 
molecule,  such  as  in  the  lone  pairs  of  a  fluorine  atom,  then  the  Fermi 
hole  may  be  more  strongly  dependent  on  the  location  of  the  probe 
electron  than  where  the  electrons  are  strongly  localized.   In  such 
cases,  the  localized  orbitals  determined  by  the  Fermi  hole  method  may 
reflect  the  locations  of  the  probe  electron  points  more  strongly  than 
they  are  reflected  in  well  localized  regions.   In  such  regions,  .however, 
there  may  be  no  physically  meaningful  way  to. distinguish  between  the 
localized  orbitals  determined  by  this  method  and  those  determined  by  any 
other  method.   In  these  situations,  the  Fermi  hole  method  may  provide  a 
practical  method  for  avoiding  the  convergence  problems  which  may  be 
expected  for  iterative  methods  when  the  electrons  are  not  well 
localized. 

The  electronic  structure  of  most  common  stable  molecules  may  be  • 
described  by  an  obvious  set  of  chemical  bonds,  lone  pair  orbitals,  and 
innershell  atomic  orbitals.   This  is  reflected  in  the  success  of  methods 
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such  as  molecular  mechanics  '   for  predicting  the  geometries  of 
complex  molecules.   Tlie  ]ocalized  orbitals  of  such  molecules  are 
unlikely  to  be  the  objects  of  much  interest  in  themselves,  but  they  may 
be  useful  in  the  calculation  of  other  properties  of  a  molecule,  such  as 

the  correlation  energy,   spectroscopic  constants,   '   and  other 

46  47 
properties.   '    The  selection  of  a  set  of  probe  electron  positions  for 

one  of  these  molecules  is  simple  and  unambiguous,  and  the  method 

presented  here  has  significant  practical  advantages  compared  to 

alternative  methods  for  transforming  canonical  SCF  molecular  orbitals 

into  localized  molecular  orbitals. 

For  some  molecules,  the  pattern  of  bonding  may  not  be  unique  or  it 

may  not  be  entirely  obvious,  even  when  the  geometry  is  known.   For 

example,  two  or  more  alternative  (resonance)  structures  may  be  involved 

in  the  electronic  structure  of  such  molecules.   The  localizpd  orbitals 

of  such  molecules  may  be  of  interest  in  themselves,  in  order  to 

characterize  the  electronic  structure  of  such  molecules,  in  addition  to 

their  utility  in  subsequent  calculations.       In  order  to  apply  the 

current  method  to  such  molecules,  the  Fermi  hole  mobility  function"*"^ '^^ 

must  be  used  to  resolve  any  ambiguities  which  may  arise  in  the  selection 

of  the  probe  electron  positions.   If  two  or  more  bonding  schemes  are 

possible,  the  positions  of  the  probe  electrons  should  be  cJiosen  to 

provide  the  minimum  values  of  the  Fermi  hole  mobility  functions  F(r)  or 

the  mobility  function  difference  F(r)-F„(p). 

In  the  case  of  methylacetylene,  for  example,  the  C-C  single  bond  may 

be  determined  by  a  single  point  midway  between  the  carbon  atoms,  where 

F(r)  is  less  than  Fq(p).   Any  attempt  to  represent  this  portion  of  the 

molecule  with  a  double  bond  would  require  placing  a  probe  electron  away 

from  the  C-C  axis,  in  a  region  where  F(r)  is  greater  than  Fp^(p). 
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Consequently,  it  is  not  possible  to  represent  methylacetylene  with  a 
structure  like  H-C=C=CH2  without  placing  one  or  more  probe  electron 
points  in  regions  where  the  Fermi  hole  is  unstable. 

In  extreme  cases,  even  the  Fermi  hole  mobility  function  may  fail  to 
provide  unambiguous  positions  for  the  probe  electrons.   This  is  expected 
in  highly  conjugated  aromatic  molecules,  metallic  conductors,  and  other 
highly  delocalized  systems.   For  these  electrons,  the  method  presented 
here,  as  well  as  all  other  methods  for  calculating  localized  orbitals, 
are  entirely  arbitrary.  The  electronic  structure  of  such  a  delocalized 
system  may  be  represented  by  an  unlimited  number  of  localized 
descriptions,  each  of  which  is  equally  valid. 

If  there  is  a  need  for  imposing  a  localized  description  on  a  highly 
delocalized  system,  the  current  method  would  be  no  less  arbitrary  than 
existing  alternatives.   The  arbitrariness  of  the  current  method  would  be 
manifested  in  the  choice  of  the  probe  electron  positions  for  the 
delocalized  electrons.  However,  the  current  method  would  continue  to 
provide  practical  advantages  over  alternative  transformations.   These 
advantages  include  the  absence  of  integrals  to  evaluate,  the  absence  of 
iteratively  repeated  calculations,  and  the  absence  of  convergence 
problems. 

Boys  Localization 

The  most  widely  used  form  of  localized  orbitals  are  those  orbitals 

5  6 
based  on  the  method  of  Boys.  '   The  integral  transformation  procedure 

in  any  localization  procedure  can  be  a  time  limiting  step.   In  the 

Edmiston-Reudenberg  method  the  two  electron  repulsion  integrals  must  be 

transformed  an  N  computational  step,  but  the  Boys  method  may  be 
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formulated  in  terms  of  products  of  molecular  dipole  integrals.   The 

dipole  integrals  are  one-electron  integrals  therefore  the  transformation 

3 
is  on  the  order  of  N  ,  making  the  Boys  localization  the  method  of 

choice.   One  disadvantage  of  this  method  is  that  it  is  an  iterative 

method  which  may  be  prone  to  convergence  difficulties. 

Methods  of  Boys  Localization 

There  exists  a  unitaiy  transformation  relating  a  delocalized  set  of 
orbitals  to  a  localized  set,  but  the  form  of  this  transformation  is  in 
general  unknown.   The  common  method  foi  solving  this  problem  is  to  do  a 
series  of  unitary  transformations  that  increase  the  degree  of 
localization  of  a  set  of  orbitals. 

Given  a  set  of  orbitals  f  one  can  increase  the  localization  by 
doing  a  unitary  transformation  U  , 

fI+1   „I.I 

f    =  U  f  (1-19) 

where  f    is  the  resulting  set  of  more  localized  orbitals.   The 
original  U  matrix  formulation  is  based  on  a  sequence  of  pairwise 
rotations,  as  proposed  by  Edmiston  and  Reudenberg.^  In  this  procedure, 
N  orbitals  are  localized  by  rotating  a  pair  of  orbitals,  then  a  second 
pair  of  orbitals  is  rotated,  .  .  .,  etc.  until  all  N(N-l)/2  pairs  of 
orbitals  have  been  rotated.   Since  the  rotations  are  done  in  a  specific 
order,  the  localized  orbitals  obtained  will  be  dependent  on  the  order  of 
rotations. 

Leonard  and  Luken  have  developed  a  second  order  method  that  does  all 
of  the  N(N-l)/2  rotations  at  once  rather  than  one  at  a  time."^^  The  use 
of  a  second  onder  method  may  have  the  additional  benefit  of  improving 
the  convergence  difficulties  encountered  in  iterative  methods.   Their 
method  is  outlined  below. 
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The  NxN  unitary  transformation  matrix  for  the  localization  can  be 
written  as 

U  =  WR.  (1-20) 

where  W  is  a  positive  definite  matrix  defined  by 

W  =  (RR^)-l/2^  (l_2i) 

The  matrix  R  will  be  defined  as 

R  =  NT,  (1_22) 

where 

T  =  1  +  t   ,  (1-23) 

and  t  is  an  antisymmetric  matrix 


t^  =  -t   . 


(1-24) 


The  N  matrix  is  a  diagonal  matrix  which  normalizes  the  columns  of  R.   By 

application  of  U  to  a  set  of  orbitals  {f^,  ,  .  .  ,  f  )  one  produces  a 

set  of  more  localized  orbitals  {f^,  .  .  .  ,  f^} .   The  new  value  of  the 

localization,  G',  is  given  by 

N 
G'  =  E  (i'i',i'i')  (1_25) 

i'=l 

The  new  f  orbitals  can  be  thought  of  in  terms  of  a  pertubative 
expansion, 

f  =  f  +  Z  t  f  +  Z  t  t  f  (1_26) 

I       I,J^  -^ 

to  second  order.   The  t^  matrix  does  2x2  rotations  that  mix  in  portions 
of  all  the  occupied  orbitals  into  orbital  f.   The  third  term  t-t.  is 
the  product  of  a  pair  of  2x2  rotations.   When  this  form  of  the  f 
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orbitals  is  substituted  into  the  G'  equation,  you  obtain 

G'  =  Gq  +  r  tjG^(I)  +  Z   tjtjG2(I,J)    .  (1-27) 

Leonard  and  Luken   include  the  G^  second  order  term  to  accelerate 
convergence  when  one  is  in  the  quadratic  region;  only  the  first  order 
term  G^  is  calculated  initially,  and  until  the  quadratic  region  is 
encountered.  In  practice,  standard  procedures  for  localization  often 
can  take  several  hundred  iterations  to  converge.   These  second  order 
procedure  described  above  seldom  takes  more  than  20  cycles,  and  the 
energy  of  orbitals  related  by  symmetry  (C-H  bonds  in  benzene,  etc.)  is 
usually  reproduced  to  seven  significant  figures. 

Results 

One  example  of  the  Boys  method  of  Leonard  and  Luken  was  given  in  the 
Fermi  hole  method  section,  in  this  section  we  will  show  the  localized 
orbitals  foi  the  Ni(CO),  ion.   In  a  recent  experimental  paper  by  Reutt 
et  al. ,  the  photoelectron  spectroscopy  of  Ni(CO),is  reported  in  order  to 
clarify  the  nature  of  the  transition  metal  carbonyl  bond.   Since  the 
spectrum  is  interpreted  in  terms  of  localized  orbitals  on  both  the  metal 
and  the  carbonyl  groups,  a  first  step  in  any  quantum  chemical  treatment 
of  the  problem  is  to  localize  the  orbitals. 

This  system  is  somewhat  complicated  because  it  has  an  unpaired 
electron.   The  three  T„  orbitals  are  occupied  by  five  electrons  leading 
to  a  triply  degenerate  ground  state.   The  MOs  are  localized  using  the 
Leonard  and  Luken  implementation  of  the  Boys'  procedure.   The 
localization  breaks  the  orbitals  into  several  classes;  (1)  oxygen  lone 
pairs,  (2)  carbon  oxygen  T  (banana)  bonds,  (3)  nickel  carbon  bonds,  and 
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(4)  nickel  d-orbitals  of  two  types  namely  E  and  T„  type  orbitals.  The 
localization  may  also  be  done  on  the  unoccupied  orbitals;  this  separates 
the  unoccupied  orbitals  into  two  sets  (1)  nickel  carbon  antibonds  and 
(2)  carbon  oxygen  antibonds. 

In  Fig.  1-5,  an  iso-value  plot  of  one  of  the  nickel  carbon  bonds  is 
shown.  The  orbital  shows  a  large  amplitude  near  the  carbon  atom, 
indicative  of  large  p-orbital  contributions  on  the  carbon  and  a 
relatively  small  contribution  from  the  nickel  d-orbitals.  As  one  can  see 
from  Fig.  1-6,  a  sizable  contribution  to  bonding  comes  from  one  of  the 
partially  occupied  nickel  d-orbitals.  The  nickel  carbon  antibond  shown 
in  Fig.  1-7  possesses  a  large  node  along  the  internuclear  axis.   The 
nickel  carbon  bond  is  expected  to  be  quite  weak,  because  it  is  composed 
of  a  sum  of  the  two  bonding  orbitals  shown  in  Figs.  1-5  and  1-6.  The 
diffuse  nature  of  the  photoelectron  spectra  indicates  the  population  of 
additional  vibrational  modes,  resulting  from  a  distortion  from  a 
tetrahedral  geometry.   The  weak  nickel  carbon  bonds  would  allow  for  such 
a  distortion  to  take  place  in  the  ion.   A  further  application  of  the 
Boys  method  would  be  to  include  the  localized  orbitals  into  a  limited  CI 
calculation  to  see  if  one  could  predict  the  photoelectron  spectra  for 
Ni(CO)^. 

The  use  of  any  localized  orbital  technique  does  not  add  or  subtract 
information  from  the  overall  wavefunction.  These  methods  only  divide 
orbitals  into  more  chemical  pieces  allowing  for  easier  interpretation  of 
experimental  results. 
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Figure  1-5:  An  iso-value  localized  orbital  plot  of  a  nickel  carbon 

bonding  orbital  in  the  Ni(CO)^  molecule.  The  dashed  lines 
indicate  an  orbital  amplitude  of  -0.05  a.u.  per  cubic  bohr. 
The  solid  lines  indicate  an  orbital  amplitude  of  0.05  a.u. 
per  bohr. 
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Figure  1-6:  An  iso-value  localized  orbital  plot  of  a  nickel  carbon  non- 
bonding  orbital  in  the  Ni(CO)^  molecule.   The  dashed  lines 
indicate  an  orbital  amplitude  of  -0.05  a.u.  per  cubic  bohr, 
The  solid  lines  indicate  an  orbital  amplitude  of  0.05  a.u. 
per  bohr. 
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Figure  1-7:  An  iso-value  localized  orbital  plot  of  a  nickel  carbon  anti- 
bonding  orbital  in  the  Ni(CO)^  molecule.   The  dashed  lines 
indicate  an  orbital  amplitude  of  -0.05  a.u.  per  cubic  bohr. 
The  solid  lines  indicate  an  orbital  amplitude  of  0.05  a.u 
per  bohr. 


CHAPTER  TWO 
LANTHANIDE  CHEMISTRY 

Background 

The  past  decade  has  seen  a  dramatic  increase  in  interest  and 
activity  in  lanthanide  and  actinide  chemistry.   Not  only  has 
considerable  knowledge  been  gained  in  the  traditional  area  of  inorganic 
/-element  chemistry,  but  much  modern  work  is  concerned  with  organo-/- 
element  reactions,   and  the  use  of  lanthanides  and  actinides  as  very 
specific  catalysts. ^°'^^  Unlike  the  corresponding  chemistry  involving 
the  d  metals,  very  little  explanation  is  offered  for  much  of  this 
chemistry. 

The  electronic  structure  of  these  systems  is  difficult  to  calculate 
from  quantum  chemical  means  for  several  reasons.  Most  of  the  complexes 
of  real  experimental  interest  are  large.   In  addition,  veiy  little  about 
/-orbitals  as  valence  orbitals  is  known,  although  experience  is  now 
being  gained  on  the  use  of  f   orbitals  as  polarization  orbitals. 
Finally,  the  /-orbital  elements  are  sufficiently  heavy  that  relativistic 
effects  become  important.   Very  few  ab  initio  molecular  orbital  studies 
have  been  reported  on  /-orbital  systems. ^^  Extended  Huckel  calculation, 
however,  have  been  successful  in  explaining  some  of  this  chemistry. ^^ 
Scattered  wave  and  DVM  Xa  studies  of  /-orbital  systems  have  also  proven 
effective,  especially  in  examining  the  photoelectron  spectroscopy  of 
reasonably  complex  systems. 

We  examine  an  Intermediate  Neglect  of  Differential  Overlaps  (INDO) 
technique  for  use  in  calculating  properties  of  /-orbital  complexes.   At 
the  Self-Consistent  Field  (SCF)  level  this  technique  executes  as  rapidly 
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on  a  computer  as  does  the  Extended  Huc.kel  method,  and  considerably  more 
rapid  than  the  scattered  wave  Xa  method.   Since  the  electrostatics  of 
the  INDO  method  are  realistically  represented,  molecular  geometries  can 
be  obtained  using  gradient  methods.^''  Since  the  INDO  method  we  examine 
contains  all  one-center  two-electron  terms  it  is  also  capable  of 
yielding  the  energies  of  various  spin  states  in  thf ^e  systems.   With 
configuration  interaction  (CI)  this  model  should  also  be  useful  in 
examining  the  UV-visible  spectra  of  /-orbital  complexes.   Preliminary 
studies  of  /-orbital  chemistry  using  an  INDO  model  have  been  disclosed 

by  Clack  and  Warren   and,  more  recently,  by  Li-Min,  Jing-Quing,  Guang- 

59 
Xian  and  Xiu  Zhen.    The  method  we  examine  will  differ  from  their 

methodology  in  several  areas,  as  discussed  below. 

Several  problems  unique  to  an  INDO  treatment  of  these  systems  must 

be  considered,  and  we  have  very  little  ab  initio  work  to  guide  us.  As 

mentioned,  what  role  do  relativistic  eflects  play?  Although  we  might 

hope  to  parameterize  scalar  contributions  through  the  choice  of  orbitals 

and  pseudo-potential  parameters,  spin  orbit  coupling,  often  larger  than 

crystal  field  effect.^,  will  need  to  be  considered  at  some  later  stage. 

Since  /-orbitals  a.  c  generally  tight,  and  ligand  field  splittings  thus 

small,  a  great  many  states  differing  only  in  their  /-orbital  populations 

lie  very  close  in  energy.  These  near  degeneracies  often  prevent 

"automatic"  SCF  convergence,  a  problem  with  which  we  must  deal  for  an 

effective  model.   The  nature  of  the  valence  basis  set  itself  is  in 

question.   Are  the  filled  5p  and  the  vacant  6p  of  the  lanthanides  both 

required  for  a  proper  description  of  their  compounds? 
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Model 


The  INDO  model  Hamiltonian  that  we  use  was  first  disclosed  by  Pople 
and  collaborators, ^°  and  then  adjusted  for  spectroscopy^^  and  extended 
to  the  transition  metal  series."-^"*  The  details  of  this  model  are 
published  elsewhere.   ~  '  To  extend  this  model  Hamiltonian  to  the  /- 
orbital  systems  we  need  first  a  basis  set  that  characterizes  the  valence 
atomic  orbitals,  and  that  is   subsequently  used  for  calculating  the 
overlap  and  the  one-  and  two-center  two-electron  integrals. 

Subsequent  atomic  parameters  that  enter  the  model  are  the  valence 
state  ionization  potentials  used  for  calculating  ont^-center  one-electron 
"core"  integrals  and  the  Slater-Condon  f""  and  G^   integrals  that  are  used 
for  thr  formation  of  one -center  two-electron  integrals.   The  evaluation 
of  these  integrals  using  experimental  information  has  traditionally  made 
this  model  highly  successful  in  predicting  optical  properties. ^^'^^ 

We  -mploy  in  this  model  one  set  of  pure  parameters,  the  resonance  or 
B(k)  parameters;  for  each  lanthanide  atom  we  decided  to  use  B(s)  =  B(p), 
B(d)  and  B(/).   These  parameters  will  be  chosen  to  give  satisfactory 
geometries  of  model  systems.   Another  choice  is  one  that  gives  good 
prediction:'  of  UV-visible  spectroscopy . "^' ^^  These  values  seldom  differ 
much  from  those  chosen  to  reproduce  molecular  geometry. 

In  this  initial  work  all  two-center  two-electron  integrals  required 
for  the  INDO  modol  Hamiltonian  are  calculated  over  the  chosen  basis  set, 
as  are  the  one-center  two-electron  V°   integrals.   An  alternate  choice 
would  be  one  that  focuses  on  molecular  spectroscopy.   In  such  a  case, 
and  one  that  we  have  to  investigate  subsequently,  the  one-center  two- 
electron  F°   could  be  chosen  from  the  Pariser  approximation^''  F°(n)  = 
IP(n)  -  EA(n),  (IP  =  Ionization  Potential,  EA  =  Electron  Affinity)  and 
the  two-electron  two-center  integral  from  one  of  the  more  successful 
function?  established  for  this  purpose. ^^"'^'^ 
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At  the  SCF  level,  we  seek  solutions  to  the  pseudo-eigenvalue  probl 


em 


F  C  =  C  £ 


(2-1) 


with  F,  the  Fock  or  energy  matrix,  C,  the  matrix  compound  of  Molecular 
Orbital  (MO)  coefficients,  and  e,  a  diagonal  matrix  of  MO  eigenvalues. 
The  above  equation  is  for  the  closed  shell  case  (all  electrons  paired). 
The  unrest)  icted  Hartree  Fock  case  is  discussed  in  detail  elsewhere, ^^ 
as  is  the  open  .hell  restricted  case.'^^"'^^  Although  nearly  all  f- 
orbital  systems  are  open  shell,  consideration  of  the  closed  shell  case 
demonstrates  the  required  theory  and  is  considerably  simpler. 
Within  the  INDO  model,  elements  of  F  are  given  by: 
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(2-2c) 


where 

(Uv|ffX)  =  JdT(l)dT(2)  )(^(1)  x^(l)  rll   x*(2)  X;^(2)  (2-3) 

P  is  the  first  order  density  matrix,  and  since  one  assumes  that  the 
Atomic  Orbital  (AO)  basis  [x^}  is  orthonormal  it  is  identical  to  the 
charge  and  bond  order  matrix,  given  by 
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with  n^  the  occupation  of  MO  #^,  n^  =  0,1,2.   In  Eq.  (2-2),  F^^  refers 
to  a  matrix  element  <xJ|f|x^>  with  AG  X^  centered  on  atom  "A".   The  core 


integral 


C  -  <^l  -  I '^  -  V  v*l<) 


(2-5) 


is  essentially  an  atomic  term  and  will  be  estimated  from  spectroscopic 

data  as  described  below.   V  is  an  effective  potential  that  keeps  the 

valence  orbital  x^  orthogonal  and  non-interacting  from  the  neglected 

inner-shell  orbitals.  The  choice  of  an  empirical  procedure  for  U^^  will 

remove  the  necessity  for  explicit  consideration  of  this  term.   The  bar 

over  an  orbital  in  an  integral, 

_A_A 
such  as  (u  y|  indicates  that  the  orbital  x^  is  to  be  replaced  with  an  s 

symmetry  orbital  of  the  same  quantum  number  and  exponent.  The 
appearance  of  such  orbitals  in  the  theory  is  required  for  rotational 
symmetry  and  compensates  for  not  including  other  two  center  integrals  of 
the  NDDO  type;^"^  i.e.  (u\^\,    X^^x^.      The  last  term  in  Eq.  (2-2a) 
represents  the  attraction  between  an  electron  in  distribution  X^*  xf  and 
all  .luciei  but  A.   The  rationale  for  replacing  integral 

(/\R-^'\A   <-—  (aVisV)  (2-6) 

is  given  elsewhere,  and  compensates  for  neglected  two  center  inner 
shell-valence  shell  repulsion^'' '^^  and  neglected  valence  orbital 
(symmetrical)  orthogoHalization. ^^^  29  ^^  .^  ^^^  ^^^^  ^^^^^^  ^^  ^^^^  ^ 
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and  is  equal  to  the  number  of  electrons  of  neutral  atom  B  that  are 
explicitly  considered;  i.e.  4  for  carbon,  8  for  iron,  4  for  cerium,  etc, 

S  of  Eq.  (2-2c)  is  related  to  the  overlap  matrix  D,"^^"-^^  and  is 
given  by 


^v  =  2  ^u(l)v(l)  ^u(l)v(l)  C^CDhd)) 
1=0 


(2-7) 


where  gy(j^)^(2)  is  the  Eulerian  transformation  factor  required  to  rotate 
from  the  local  diatomic  system  to  the  molecular  system,  (y(l)|v(l))  are 
the  Sigma  (1=0),  pi(l=l),  delta(l=2)  or  phi  (1=3)  components  to  the 
overlap  in  the  local  system,  and  /^(j^^^)  are  empirical  weighting 
factors  chosen  to  best  reproduce  the  molecular  orbital  energy  spread  for 
model  ab  initio  calculations.   We  have  made  little  use  of  this  /  factor, 
and  set  all  /  =  1  except  between  p  symmetry  orbitals  viz.^-^"^"^'^^ 


"  ""-^^^hnvn   ^P^  IP^  ) 


(2-8) 


Basis  Set 


In  general  ZDO  methods  choose  a  basis  set  of  Slater  Type  Orbitals 


(STO) 


nlm 


12C1 


2n+l 


2n! 


1/2 


r"-V^^  Yj  (9,*) 


(2-9a) 


where  X^{Q,^)   are  the  real,  normalized  spherical  harmonics.   Atomic 

orbitals  X  are  expressed  as  fixed  contractions  of  these  fR   ] 
^  nlm-" 
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^i  "  ^  ^nlm  ^nlm  (2-9b) 

In  general  a  single  R^^^  function  describes  the  s  and  p  orbitals  for 
most  atoms.   The  d  orbitals  of  the  transition  metals,  however,  require 
at  least  a  double-C  type  function  (two  terms  in  2-9b)  for  an  accurate 
description  of  both  their  inner  and  outer  regions.   For  the  lanthanides 
wo  have  examine',1  basis  sets  suggested  by  Li  Le-Min  et  al.,^^  by  Bender 
and  Davidson,   and  by  Clementi  and  Roetti.^^  In  the  latter  case,  the 
two  major  contributors  of  Eq.  (2-9b)  in  the  valence  orbitals  of  the 
double-C  atomic  calculations  were  selected,  and  these  functions  were 
renormalized  with  fixed  ratio  to  yield  the  required  nodeless  double-C 
functions  for  illDO.   We  were  unable  with  any  of  these  choices  to  develop 
a  sy-ctematic  model  useful  for  predicting  molecular  geometries  (see  later 
discussion  of  resonance  integrals). 

We  have  adapted  the  following  procedure  on  selecting  an  effective 

1   .      80  " 

basis  set.    Knappe  and  Rosch  calculated  the  lanthanides  and  their 

mono-positive  ions  using  the  numerical  Dirac-Fock  relativistic  atomic 

81 
program  of  Desclaux.    From  these  wavefunctions  j ndial  expectation 

2        "^ 
values  <r>,  <r  >  and  <r  >  are  calculated  for  6s,  6p,  5d  and  4/ 

functions.   The  6s,  5d  and  4/  wavefunctions  were  obtained  by  Dirac-Fock 

calculations  on  the  promoted,  Af-^5d^6s^   configuration;  the  6p  from 

calculations  in  which  a  5d  electron  was  promoted,  ^f~^6s^6p^. 

Wavefunctions  for  the  mono-positive  ions  are  obtained  from  4/'~^3d^6s-^ 

and  4/"  6-:  6p  respectively.   A  generalized  Newton  procedure  was  then 

used  to  determin  exponents  (C)  and  coefficients  a^^^  for  a  given  set  of 

<r>,  <r  >  and  <r  >  with  functions  of  the  form  of  Eq.  (2-9b).   Again,  as 

in  the  transition  metal  atoms,  we  found  that  a  single  C  function  fits 

the  ns  and  np  atomic  functions  well  in  the  regions  where  bonding  is 
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important,  but  the  (n-l)d,  and  now  the  (n-2)/  require  at  least  two  terms 
in  the  expansion  of  Eq.  (2-9b).   This  is  demonstrated  for  the  Ce^  ion  in 
Figure  2-1,  where  it  is  shown  that  a  single-^  expansion  is  poor  for  the 
outer  region  of  the  4/  function. 

In  Figure  2-2  the  value  of  <r>  is  plotted  versus  atomic  number.   The 
contraction  of  the  6s  and  6p  orbitals  due  to  relativistic  effects  (DF 
vs.  HF)  is  quite  apparent  here,  and  is  a  consequence  of  the  the  greater 
core  penetration  of  these  orbitals.   Subsequent  expansion  of  the  4/  and 
5d,  now  with  increased  shielding,  results.   After  some  experimentation 
we  use  the  Dirac-Fock  values  obtained  from  th.:  mono-positive  ions.   The 
basis  set  adopted  is  given  in  Table  2-1. 

The  4/  and  5d  functions  are  quite  compact.   At  typical  bonding 
distance  (4/4/|uu)  and  (5d5d  |uu)  are  essentially  R'^.  Because  of  this 
we  calculate  all  two-cenrer  two-electron  integrals  with  the  ^  values  in 
Table  2-1.   This  value  is  chosen  to  match  the  accurate  F°  Slater-Condon 
Factors  obtained  from  the  numtirical  atomic  calculations  by  a  single 
exponent,  via 

F°(4/4/)  =  0.200905  C  (4/)  (2-10a) 

F^CSdSd)  =  0.164761  C  (5d)  (2-lOb) 

F-'CesGs)  =  0.139803  C  (6s)  (2-lOc) 

F°(6p6p)  =  0.139803  C  (6p)  (2-lOd) 

The  error  in  ca]culating  two-centered  two-rlectron  integrals  at  typical 

bonding  distances  with  this  single-C  approximation  is  well  under  V/,,    and 

this  procedure  is  much  simpler. 

Core  Integrals 

The  average  energy  of  a  configuration  of  an  atom  or  ion  is  given 
by«2,83 
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Single  vs.  Double  Zeta  4f-ST0  Orbital  Amplitude 
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Figure  2-1:   Radial  wavefunction  for   the  4/  orbital  of  Ce     with  single-C 
and  double-C  Slater   type  orbitals   (STOs). 
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DF  vs.  HF  Average  Valuas  of  r  for  th«  4f,5d,6»  and  6p  Orbitoia 
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Figure  2-2:  Average  value  of  r  for  the  valence  orbitals  of  the 

lanthanides  from  a  relativistic  calculation  (DF)  and  a  non- 
relativistic  calculation  (HF). 
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with  W^j,  the  average  two  electron  energy  of  a  pair  of  electrons  in 
orbitals  X^  and  X-  given  by 

"ss  =  '^°(^^>  (2-12) 

Wpp  =  F'CPP)  -  2/25f2(pp) 

W^j  =  F°(dd)  -  2/63F^(dd)  -  2/63  F^(dd) 

W^^  =  F^iff)   -   4/]95F^(/J=)  -  2/143  F^(ff)   -   100/5577  F^(//) 

^sp  =  ^°^^P^  -  1/6g\sp) 

^sd  =  ^°^^^^  -  ^/^°  ^^(sd) 
Wg^  =  F'Csf)  -  1/lA  G^(s/) 

V  "  ^°^P^^  -  1/1^  G^(pd)  -  3/70  G^(pd) 
Wp^  =  F°(pj=)  -  3/70  G^p/)  -  2/63  G^(pj=) 
W^j.  =  F°(dj:)  -  3/70  G^(d/)  -  2/105  G^(d/)  -  5/231  G^(d/) 
The  core  integrals  U.^,  Eq.  (2-5),  are  then  evaluated  by  removing  an 
electron  from  orbital  X • ,  and  equating  the  difference  in  configuration 
energy  between  cation  and  neutral  to  the  appropriate  observed  IP(n).   We 
prefer  this  procedure  rather  than  that  suggested  by  others  that  average 
the  value  obtained  from  IP(n)  and  EA(n) . ''^'^^'^'^ 

There  are  a  great  many  low  lying  configurations  of  the  lanthanide 

atoms  and  their  ions.   The  lowest  terms  of  Ce,  Gd  and  Lu  come  from  f'^d 

2 
s  ,  while  the  remaining  lanthanide  atoms  have  the  structure  /""^  s^. 

Two  processes  are  then  possible  for  6s  electron  ionization: 

I  f-\h^   -  /"-^d^s^   .  (s) 

II  f-'     s'   .       f-h'      .    (s) 


53 
The  ionization  energy  of  a  6s  electron  from  I  is  systematically  0.4  - 
0.5  eV  larger  I  lian  that  obtained  from  II.   When  combined  with  Eq.  (2- 
11),  the  estimate  for  U^^  differ  by  less  rhan  0.1  eV.  That  is,  choosing 
the  values  of  process  I,  the  use  of  Eq.  (2  11)  predicts  the  values  of 
process  II  within  0.1  eV.   We  thus  choose  the  values  of  process  I  shown 

in  Table  2-2.   These  values  are  obtained  from  the  promotion  energies  of 

85   8fi 
Brewer   '    and  then  smoothed  by  a  quadratic  fit  thrMUgliout  the  series. 

For  completeness,  we  also  give  the  values  of  process  II. 

The  lowest  configuration  containing  a  5d  electron  is  j^'^d^s^ 

throughout  the  series,  and  5d  ionizations  are  obtained  from 

III  f-'  ih'   .  f-h'     .   (d) 

The  ionization  potentials  for  the  6p  can  be  obtained  from  two 
processes: 

IV  /"-2  s^p^  .  f-s'^        .  (p) 

V  f-^  s'p   .  f-V     .   (p) 

Ionization  from  process  IV  is  nearly  constant  at  3.9  eV,  from  V  at 

ji  -2  2 
4.6  eV.   The  f  ''s  configuration  is  lower  for  all  the  lanthanides 

except  Ce(/dsp),  and  Gd  and  Th{f^~^s^\>).      Using  1 1.-  ionization 

potentials  of  process  IV,  and  Eq.  (2-11),  we  predict  the  values  of 

process  V  to  within  0.2  eV.   Ue  do  not  consider  this  error  significant, 

and  thus  use  the  smoothed  values  from  IV  given  in  Table  2-2.   The  values 

from  process  V  are  also  given  in  the  table  for  comparison. 

For  a  /  orbital  ionization,  we  consider  the  two  pro'-esses 

VI  f~^ds'   ^f-W   .   (f) 

VII  /-2  s^  .  f-V     .   (/) 

(compare  with  I  and  II).   As  seen  jn   Table  2-2  the  values  form  the  t.wo 
processes  are  very  different.   From  Eq.  (2-11) 

U^^  (VI)  =  IP(VI)  -  (m-4)  V^^  -  2W^^  -  W^^  (2- 13a) 

U^^  (VII)  =  IP(VII)  -  (m-3)  W^^  -  2W^^  (2-13b) 
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Unlike  the  analogous  situation  for  the  6s  and  6p  orbitals,  use  of 
Eq.  (2-13a)  to  find  U^^,  and  use  of  this  value  in  Eq.  (2-13b)  to  predict 
IP(VII)  is  not  successful,  and  would  require  the  scaling  of  the  large 
F°(/f)  integral  often  performed  in  methods  parameterized  on  molecular 
spectroscopy. ^^'^-^'^^ 

As  with  the  transition  metal  nd  orbitals  ve  might  envision  the 
following  procedure.  We  assume  that  the  lanthanide  atom  in  a  molecule 
is  a  weakly  perturbed  atom.   The  lowest  energy  configuration  of  the  atom 
should  than  be  most  important  in  determining  U„,.   Ue  create  a  two-by- 
two  interaction  matrix 


fE(f-W)   -   X     ^ 


Eif-hh  -  xj 


l^2J 


0 


(2-14a) 


where  V  is  an  empirical  mixing  parameter,  and  C^^  and  C^^  determines  the 
relative  amounts  of  each  of  the  two  configurations  that  are  important. 
The  exact  value  of  V  would  depend  on  a  given  molecular  situation.   C  ^ 
is  then  given  by 


''    '  ux' 


(2-14b) 


-3^  2, 


X  =  C  /C  -  E(f""'d  s-^)  -  E(f"-V^ 


(2-14c) 


(E(f-h  sh   -  Erf"-^.^ 


n2 


2V 


The  values  of  C^  appear  in  Table  2-2,  where  we  have  used  the  values  of 
E(f~   ds  )  and  ECJ^'^s^)  obtained  for  the 


le  promotion  energies  of 
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Brewer   and  a  fixed  value  of  V  =  0.02  au.  Then  U^^  could  be  obtained 
from 
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"//  =  ^1^  V^'^  '   4^  V''^  (2-15) 

In  the  case  of  the  3d  orbitals  this  valence  bond  mixing  between  3d""^s^ 

and  3d   s  was  important  in  obtaining  reasonable  geometric 

63 
predictions,   an  observation  now  confirmed  in  careful  ab  initio 

87 
studies.    For  the  lanthanide  complexes  of  this  study  the  Af   orbitals 

are  quite  compact,  and  this  valence  bond  mixing  does  not  greatly  affect 
geometries.   However,  the  calculation  of  ionization  potentials  that 
result  in  states  with  reduced  /-orbitals  occupation  is  influenced. 
There  are  many  refinements  one  can  make  in  the  formation  of  a 
"mixing"  matrix  such  as  Eq.  (2-Ua).   One  might  be  to  make  V  dependent 
on  the  calculated  population  of  the  4/  and  5d  atomic  orbitals.  However, 
the  values  of  the  promotion  energies  we  obtain  from  Brewer  are  so 
different  than  those  that  we  obtain  from  oui  own  numerical  calculations 
on  the  average  energy  of  a  configuration,  Table  2-3,  that  for  the  moment 
we  choose  a  76%  :  2A%   mix  of  E(f-^ds^)    :    Eif'^s^)    for  all  the  atoms  of 
the  series.   This  mix  gives  reasonable  geometries  and  ionization 
potentials  for  all  molecules  of  this  study.   Further  refinements  will 
require  more  accurate  atomic  promotion  energies  and  numerical  experience 
with  the  model. 
Resonance  Parameters,  B(k) 

Each  lanthanide  atom  has  three  B  (k)  values,  B(s)  =  B(p),  B(d)  and 
B(/),  and  those  we  choose  are  summarized  in  Table  2-4.   They  are 
obtained  by  fitting  the  geometries  of  the  trihalides,  and  the  more 
covalent  bis-cyclopentadienyls  to  be  reported  elsewhere. 

Bond  lengths  are  most  sensitive  to  B(d)  and  bond  angles  to  B(p). 
These  angles  can  be  r-prcduced  solely  on  a  basis  set  including  6p 
orbitals,  and  we  have  been  able  to  obtain  satisfactory  comparisons  with 
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Table  2-3  :  Average  configuration  energy  from  Dirac-Fock  calculations^ 
on  the  J   ds  and  the  /"  s  configurations  for  all  the 
lanthanide  atoms. 


Atom 

Average 
Configuration 
Energy 

f-Vs' 

f-2^i 

Ce 

-8853.71494569 

-8853.64980000 

Pr 

-9230.41690970 

-y.;30. 37981848 

Nd 

-9616.94751056 

-9616.93446923 

Pm 

-10013.4526061 

-10013.4606378 

Sm 

-10420.0710475 

-10420.0976615 

Eu 

-10836.9533112 

-10836.8834715 

Gd 

-11264.0945266 

-11264.0439334 

Tb 

-11701.7877496 

-11701.7482691 

Dy 

-12150.1565528 

-12150.1286785 

Ho 

-12609.3663468 

-12609.3484161 

Er 

-13079.5686394 

-13079.5585245 

Tm 

-13560.9236801 

-13560.9201649 

Yb 

-14053.5770354 

-14053.5786047 

Lu 

-14557.7153258 

a)  Reference  81. 
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Table  2-4  :Resonance  integrals  (B  values)  for  the  Lanthanide  atoms  in  e 
V.  The  beta  for  the  s-orbital  is  set  equal  to  the  beta  for 
the  p-orbital. 


Atom 

B(s) 

B(p) 

B(d) 

B(/) 

Ce 

-8.00 

-8.00 

-17.50 

-80.00 

Pr 

-7.61 

-7.61 

-17.58 

-80.00 

Nd 

-7.23 

-7.23 

-17.65 

-80.00 

Pm 

-6.85 

-6.58 

-17.73 

-80.00 

Sn 

-6.46 

-6.46 

-17.81 

-80.00 

Eu 

-6.08 

-6.08 

-17.88 

-80.00 

Gd 

-5.69 

-5.69 

-17.96 

-80.00 

Tb 

-5.31 

-5.31 

-18.04 

-80.00 

Dy 

-4.92 

-4.92 

-18.11 

-80.00 

Ho 

-4.54 

-4.54 

-18.19 

-80.00 

Er 

-4.15 

-4.15 

-18.27 

-80.00 

Tm 

-3.77 

-3.77 

-18.35 

-80.00 

Yb 

-3.38 

-3.38 

-18.42 

-80.00 

Lu 

-3.00 

-3.00 

-18.50 

-80.00 
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experiment  without  the  necessity  of  including  the  5p  orbitals.  On  the 
other  hand,  orbitals  of  p  symmetry  do  seem  to  be  required  for  accurate 
predictions  of  geometry.   '^^ 

It  has  been  argued  that  the  4/  orbitals  are  not  used  in  the  chemical 
bonding  of  those  complexes  except  in  the  more  covalent  cases. ^^'^^  From 
the  present  study  we  are  lead  to  the  conclusion  that  some,  albeit  small, 
contribution  is  required  of  these  orbitals  to  obtain  the  excellent 
agreement  between  experimental  and  calculated  bond  lengths  for  the 
series  MF^,  MCl^,  MBr^  and  MI^  and  for  the  comparative  values  obtained 
for  CeF^  and  CeF^.   This  is  indicated  in  Table  2-A  by  the  large  values 
of  |B(f)|.   The  latter  values  are  a  consequence  of  the  fact  that  the  /- 
orbitals  are  tighter  than  one  usually  expects  for  orbitals  important  in 
chemica]  bonding.  Use  of  5d  orbitals  alone  will  predict  the  trends  in 
these  two  series,  but  underestimates  the  range  of  values  experimentally 
observed. 
Two  Electron  Integrals 

Several  different  interpretations  have  been  given  to  the  INDO 
scheme.   The  simplest  of  those  schemes  is  to  include  only  one-centered 
integrals  of  the  Coulomb  or  exchange  type 

(uy|vv)  or  ()jv|\)vi) 
For  an  s,p  basis  these  are  complete.   For  an  s,p,d  or  s,p,d,f  basis  they 
are  not,  and  the  omission  of  the  remaining  integrals  will  lead  to 
rotational  variance.   To  restore  rotational  invariance,  integrals  of 
this  type  might  be  rotationally  averaged, ^^  but  from  a  study  of  spectra 
it  appears  that  all  one-cpnter  integrals  should  be  evaluated.  ^"^  For 

example,  in  the  metallocenes  the  integral  (d  2  2d   Id   d   ^  i^ 

'  X  -y   xyl  xy  '^yz''  ^   •■ 

required  to  separate  the  two  transitions  that  arise  from  the  e   (d)  ^ 
e2g(d)  transitions  that  lead  to  the  ^E^^  and  ^E^     excited  states.   In 
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addition,  it  appears  that  the  inclusion  of  all  one-center  integrals 

Improves  tho  predictions  of  angles  about  atoms  with  s,p,d  basis 

89  90 
sets   '   and  considerably  improves  the  predictions  of  angles  about  the 

lanthanides.   For  these  reasons  we  include  all  the  one-center  two- 
electron  integrals.   Since  the  INDO  programs  we  use  process  integrals 
and  their  labels  in  the  MOLECULE  Format^^  only  the  additional  integrals 
need  be  included.  These  integrals  are  generated  in  explicit  form  via  a 
computer  program  that  we  have  used  in  the  past^"^  and  they  have  also  been 
recently  published  by  Schulz  et  al.^^  To  our  knowledge  all  these 
integrals  do  not  appear  in  the  literature  for  s,p,d  and  /  basis, 
although  we  have  checked  those  of  (uu|v\;)  and  (uv|\)u^  against  the 
formulas  of  Fanning  and  Fi tzpatrick.^^ 

Integrals  of  the  form  (uw|vv)  and  (uvjvu)  can  be  obtained  thi uugh 
atomic  spectroscopy,  and  their  components,  F^   and  G^,  evaluated  via 
least  square  fits 

(yy|vv)  =  E  \  F*" 

k 
(Uv|vu)  =2^2*^ 

k 

k      k 
These  F  and  G  can  then  be  used  to  evaluate  all  integrals  of  the  "F"  or 

"G"  type,  even  those  that  do  not  appe.u"  in  atomic  spectra  because  of 
high  symmetry  (i.e.  ('i^2_y2  d^Jd^^  d^^)).  Integrals  of  the  "R"  type, 
however,  cannot  be  evaluated  in  this  mannei ;  viz.(sd|dd),  (sp|pd), 
(sdipp),  (sd|//),  (s/|d/),  (ppIp/),  (dd|p/),  (prl|d/),  (sd|p/),  (pd|s/), 
(sp|d/),  and  (pf\ff).      For  this  reason  we  evaluate  all  one-center  two- 
electron  integral-  of  the  lanthanides  using  the  basis  set  of  Table  2-1, 
which  yields  the  exact  F°  value  obtained  from  the  Fock-Uirac  numerical 
calculations.  All  F*^,  G^  and  R^  integral  for  k  >  0  are  then  scaled  by 
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2/3.   This  value  of  the  scaling  is  obtained  from  a  comparison  of  the 
calculated  and  empirically  olnained^^'^^"^^  F^(ff),    F^(ff)   and  F^(//) 
values  that  implies  0.66  ±  0.U4.   Empirically  obtained  values  of  G^(/d) 
and  F  (/d)  are  far  more  uncertain  and  are  much  smaller,  and  are  thus  not 
used  to  obtain  this  scaling  value  between  calculated  and  experimental 

values. 

i 

At  this  point  it  seems  appropriate  to  point  out  the  differences  of  the 
present  INDO  model  to  that  suggested  by  Li  Le-Min  et  al.^^  In  the 
latter  formalism  only  the  conventional  one-center  two-electron  integrals 
are  included  leading  to  rotational  variance.   In  addition,  the 
Wolfsberg-Helmholz  approach  is  used  for  the  resonance  integral  B, 
B^j=(IP(i)+IP(j))S^ V2.   No  geometry  optmization  has  been  reported 
within  their  model.  ^  Further  differences  are  the  restriction  to 
single-C  STOs  and  the  smoothing  of  the  valence  orbital  ionization 
potentials  for  the  lanthanides  via  Anno-type  expressions.^^ 


Procedures 


The  input  to  the  INDO  program  consists  of  molecular  coordinates  and 
atomic  numbers.  Molecular  geometries  are  obtained  automatically  via  a 
gradient  diivpn  quasi-Newton  update  procedure, ^^  using  either  the 
restricted  or  unrestricted  Hartree  Fock  formalism.  All  UHF  calculations 
are  followed  by  simple  annihilation.^^ 

Self-consistent  field  convergence  is  a  problem  with  many  of  these 
systems.   For  this  reason  electrons  are  assigned  to  molecular  orbitals 
that  are  principally  f   in  nature  according  to  the  number  of  /-electrons 
in  the  system,  and  the  symme^^y  of  the  system.   Orbitals  with  large 
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lanthanide  5d  character  are  sought  and  assigned  no  electrons.   A 
procedure  is  then  adopted  that  extrapolates  a  new  density  for  a  given 
Fock  matrix  based  on  a  Mulliken  population  analysis  of  each  SCF  cycle.  ^-^ 

Often  this  procedure  is  not  successful.   In  such  cases  all  / 
orbitals  are  considered  degenerate,  and  they  are  equally  occupied  in  the 
highest  spin  configuration  using  the  RHF  open  shell  method.''''"  These 
vectors  (orbitals)  are  then  stored,  and  the  SCF  repeated  with  the 
specific  /  orbital  assignments  as  described  above. 

In  cases  of  slow  convergence,  a  singles  or  small  singles  and 

doubles,  CI  is  performed  to  check  the  stability  of  the  SCF,  and  the 

appropriateness  of  the  forced  electron  assignment  to  obtain  the  desired 
state.96 

Results 

The  geometries  of  CeCl^  and  LuCl^  were  used  to  determine  an  optimal 
set  of  resonance  integrals  and  conf igurational  mixing  coefficients.  No 
further  fitting  was  performed,  and  thus  the  structures  of  all  other 
compounds  are  "predictions".   The  resonance  parameters  for  the  other 
lanthanides  were  determined  by  interpolation  from  the  values  for  Ce  and 
Lu  (see  Table  2-4).   The  INDO  optimized  geometries  as  well  as  the 
remaining  cerium  and  lutetium  trihalides  are  listed  in  Table  2-5.   In 
addition  to  the  trihalides  reported,  the  geometry  of  CeF  is  also  listed 
in  Table  2-5.   One  can  see  the  agreement  with  experiment  is  good  in  all 
cases. 

The  potential  energy  of  the  trihalides  as  a  function  of  the  out  of 
plane  angle  is  very  flat.   Although  we  have  optimized  all  structures 
until  the  gradients  are  below  10~^  a.u./bohr,  the  angles  are  converged 
only  to  ±3°.   We  note,  however,  that  all  are  predicted  non-planar,  in 
agreement  with  experiments. ^"^"^'^^'^^ 


J' 
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Table  2-5  :  Geometry  and  ionization  potentials  for  Cerium  and 
Lutetium  trihalides.   Cerium  tetraf luoride  is  also 
included  in  this  table.   The  bond  distances  are  given 
in  angstroms^  angles  in  degrees  and  IPs  in  eV. 
Experimental  results  are  also  shown  where  available. 


Molecule 

Bond 
Distance 

Bond 
Angle 

Ionization 
Potential 

TNDO 

Exp. 

INDO 

Exp. 

INDO 

Exp. 

CeF^ 

2.204 

2.180 

106.8 



8.4 

8.0 

CeCl^ 

2.570 

2.569 

115.6 

111.6 

10.0 

9.8 

CeBr^ 

2.668 

2.722 

115.8 

115.0 

9.6 

9.5 

Cel^ 

2. 844 

2.927 

119.8 



9.9 



CeF, 
4 

2.099 

2.040 

109.5 

109.5 





LuF^ 

2.045 

2.020 

107.4 



b 

19.0 

LuCl^ 

2.415 

2.417 

108.2 

111.5 

18.6 

(17.4  -  18.7) 

LuBr^ 

2.528 

2.561 

108.6 

114.0 

17.8 

(16.8  -  18.4) 

LUI3 

2.726 

2.771 

115.6 

114.5 

17.7 

(16.2  -  18.1) 

a)  References  54b,  97  and  98.   Estimated  values  for  CeF^,CeI.,  and 
LuF^  from  Ref  103.  3'   3' 

b)  Thf  SCF  calculation  on  the  ion  of  LuF^  would  not  converge 
therefore  no  IP  is  reported.        ■' 
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The  experimental  range  of  the  bond  lengths  from  LnF^  to  Lnl-  is 
greater  than  we  calculate.   Our  predicted  values  for  the  trif.luorides 
and  trichlorides  are  in  good  agreement,  while  bond  lengths  for  the 
tribromides  and  triiodides  are  too  short.   Since  these  are  the  more 
polarizable  atoms  it  is  possible  that  configuration  interaction  will 
have  its  largest  affect  on  these  systems.   The  calculated  change  in  bond 
length  of  0.11  A  in  going  from  CeF^  to  CeF^  is  also  smaller  than  the 
O.IA  A  observed. 

Ionization  potentials  (IPs)  are  also  reported  in  Table  2-5.  In  all 
cases  the  INDO  values  fall  within  the  experimental  ranges.   These  values 
are  calculated  usinfr  the  ASCF  method,  and  only  the  first  IP  is 
calculated.   Experimentally  ^"^'^^  these  valued  are  somewhat  uncertain, 
but  they  are  split  by  both  crystal  field  effects,  and  by  the  large  spin- 
orbit  coupling  not  yet  included  in  our  calculations.   However,  the 
latter  interaction  is  treated  implicitly  in  the  DVM  Xa  calculations^'^'^ 
based  on  the  Dirac  equation.   Therefore,  the  Xa  result  for  the 
ionization  potentials  show  better  agreement  with  the  experiment  in  this 
aspect,  but  it  is  quite  remarkable  that  the  present  INDO  approach  is 
able  to  reproduce  the  experimental  trend  in  the  first  IP  of  the  series, 
CeX^,  X  =  F,  Cl,  Br  with  a  maximum  value  for  the  chloride,  a  feature 
noticeably  missing  in  the  DVM  Xa  results. ^'^'^ 

The  initia.1  success  of  the  INDO  model  as  implemented  here  lead  us  to 
calculate  both  geometries  and  IPs  for  the  remaining  lanthanide 

trichlorides.   These  results  are  shown  in  Table  2-6.   The  experimental 

54b  97  98 
geometries   '   '   are  very  well  reproduced  by  the  INDO  calculations. 

The  INDO  IPs  reproduce  the  characteristic  "U"  pattern  of  the  lanthanide 

atoms,  and  fall  within  the  experimental  ranges. 
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Table  2-6  :  Geometries  and  Ionization  Potentials  (IPs)  for  the 

lanthanide  trichlorides.  Bond  distances  are  reported  in 
angstroms,  bond  angles  in  degrees  and  IPs  in  eV. 
Experimental  results  are  also  given  where  available. 


Atom 

Bond 
Distance 

Bond 
Angle 

Ionization 
Potential 

INDO 

Exp. 

INDO 

Exp. 

INDO 

Exp. 

Ce 

2.570 

2.569 

115.6 

111.6 

10.0 

9.8 

Pr 

2.566 

2.553 

108.5 

110.8 

11.8 

(10.9-11.2) 

Nd 

2.563 

2.545 

112.7 



13.3 

12.0 

Pm 

2.556 



112.7 



14.4 



Sm 

2. 544 



113.0 



15.3 

(13.7-17.0) 

Eu 

2.532 



113.2 



16.4 

Gd 

2.514 

2.489 

110.0 

113.0 

17.7 

(15.5-16.5) 

Tb 

2.496 

2.478 

109.8 

109.9 

13.0 

(13.0-20.5) 

Dy 

2.479 



110.1 



14.3 

(14.0-20.0) 

Ho 

2.464 

2.459 

112.0 

111.2 

15.0 

(15.5-20.0) 

Er 

2.448 



110.9 



15.6 

(11.5-16.0) 

Tm 

2.430 



108.5 



15.9 

(15.3-21.0) 

Yb 

2.421 



109.6 



15.9 

(15.5-21.0) 

Lu 

2.415 

2.417 

108.2 

111.5 

18.6 

(17.4-18.7) 

a)  References  54b,  97  and  99, 
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To  test  the  applicability  of  our  model  to  lanthanide  atoms  not 
formally  charged  +3,  we  calculated  the  geometries  and  IPs  for  SmCl„, 
EUCI2  ^"^  ^''^"'"2  ""o^scules.   The  results  are  given  in  Table  2-7.   The 
INDO  model  gives  optimized  geometries  that  are  bent  and  in  good 
agreement  with  experimental  results. -^^^  We  note  that  this  bending  is  a 
result  of  a  small  amount  of  p-orbital  hybridization.   It  is  not 
necessary  to  invoke  London  type  forces,  and  thus  correlation,   to 
explain  this  effect. 

We  chope  Ce(NO^)g  as  our  last  example  because  it  is  one  of  the  few 
known  examples  of  a  twelve  coordinate  metal.   The  optimized  geometry  is 
summarized  in  table  2-8  and  a  plot  of  tlie  optimized  geometry  is  shown  as 
figure  2-3.  As  one  can  see  from  Table  2-8  INDO  predicts  a  geometry  that 
is  in  excellent  agreement  with  the  experimental  crystal  structure.  "'•'^^ 
Table  2-9  shows  a  population  study  of  this  complex.   AltdouKh  there  is 
some  /-orbital  participation,  it  appears  that  this  unusual  twelve 
coordinate  T^^  structure  results  from  electrostatic  forces  between  the 
ligands  and  the  relatively  large  size  of  the  Ce(IV)  ion.         '    \ 
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Table  2-7  :  Geometry  and  ionization  potential  for  SmCl^,  EuCl„,  and 

YbCl2.   Bond  distances  are  given  in  angstroms,  bond  angles 
in  degrees,  and  ionization  potentials  in  eV.  Experimental 
results  are  listed  where  available. 


Molecule 

Bond 
Distance 

Bond 
Angle 

Ionization 
Potential 

INDO 

Exp. 

I  MOO 

Exp. 

INDO 

Exp. 

SmClo 

2.584 



143.3 

130+15 

5.3 



EuCI,, 

2.576 



143.2 

135±15 

6.6 



YbCl2 

2.400 



120.2 

126±05 

3.2 



a)  Reference  100. 
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Table  2-8  :   Average  bond  distances  and  bond  angles  for  Ce(m^)Z     ion 
INUO  optimized  geometry  and  the  X-ray  crystal  structure^. 
Distances  are  in  angstroms  and  angles  in  degrees.   The  c 
subscript  on  the  oxygen  atoms  denotes  the  that  oxygen  is 
bonded  to  the  cerium  and  the  n  subscript  signifies  a  non- 
bonded  oxygen. 


Geometric 
Parameter 

INDO 

Exp. 

r(Ce-O^) 

2.55A 

2.508 

r(N-O^) 

1.256 

1.282 

r(N-O^) 

1.237 

1.235 

'3(0^-N-0^) 

121.5 

114.5 

e(O^-Ce-O^) 

50.9 

50.9 

a)  Reference  101. 
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_2 
Figure  2-3:  Plot  of  the  twelve  coordinate  CeCNO^),   ion.  Nitrogens  2,  18 

and  22  are  above  the  plane  of  the  paper,  while  nitrogens  6, 

10  and  lA  lie  below  the  plane  of  the  paper. 
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Table  2-9  :  Population  analysis  of  Ce(NO-).  .  The  oxygen  atoms  that  are 
coordinated  to  the  cerium  are  indicated  by  0  .The  Wyberg 
bond  index  is  also  given.  A  Wyberg  index  of  1.00  is 
characteristic  of  a  single  bond. 


Atom 

0 

rbital 

Atomic 
Population 

Spin 
Density 

Total 
Valence 

Ce 

s 

0.20 

0.00 



P 

0.30 

0.00 



d 

1.32 

0.00 



/ 

1.10 

1.00 



Net 

1.08 

1.00 

4.80 

N 

Net 

0.59 

0.00 

3.78 

0 
c 

Net 

-0.40 

0.00 

1.58 

0 

Net 

-0.48 

0.00 

1.80 

Bond 

Wyberg 
Bond  Index 

Ce  -  0 
c 

0.40 

N  -  0 
c 

1.37 

N  -  0 

1.22 

CONCLUSIONS 

We  develop  an  Intermediate  Neglect  of  Differential  Overlap  (INDO) 
method  that  includes  the  lanthanide  elements.   This  method  uses  a  basis 
set  scaled  to  rcjproduce  Dirac-Fock  numerical  functions  on  the  lanthanide 
mono-cations,  and  is  characterized  by  the  use  of  atomic  ionization 
information  for  obtaining  the  one-center  one-e]  r-ctron  terms,  and 
including  all  of  the  two-electron  integrals.  This  latter  refinement  is 
required  for  accurate  geometric  predictions,  some  of  which  are 
represented  her'?,  and  for  accurate  spcrtroscopi'>  predictions,  to  be 
reported  latter. 

We  have  applied  this  method  to  complexes  of  the  lanthanide  elements 
with  the  halogens.   The  geometries  calculated  for  these  complexes  are  in 
good  agreement  with  experiment,  when  experimental  values  are  available. 
The  trihalides  are  calculated  to  bu  pyramidal  in  agreement  with 
observation.   The  potential  for  the  umhre]la  mode,  however,  is  very 
flat.   The  dichlorides  of  Sm,  Eu  and  Yb  are  all  predicted  to  be  bent 
even  at  the  SCF  level,  again  in  agreement  with  experiment.   This  bending 
is  caused  by  a  small  covalent  mixing  of  ungerade  6p  and  4/  orbitals,  and 
one  need  not  invoke  London  forces  to  explain  this  observation.   Again 
the  potential  for  bending  is  very  flat. 

Within  this  model,  /-orbitals  participation  in  the  bonding  of  those 
ionic  compound  through  covalent  effects  is  small.  Nevertheless  f- 
orbitals  participation  does  contribute  to  the  pyrimidal  geometry  of  .the 
trihalides  and  the  bent  stiuctuie  of  the  dihalides.   In  addition, 
although  the  trend  of  bond  lengths  within  the  series  LnF  ,  LnCl^,  LnBr 
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and  Lul^,  and  CeF^  and  CeF,  are  reproduced  without  /-orbital 
participation,  the  range  of  values  calculated  is  considerably  improved 

when  /-orbitals  are  allowed  to  par rici pate.   For  the  twelve  coordinate 

-2       102 
CeCNO^)^  complex    reported  here,  /-orbital  participation  appears 

minor.   A  stable  complex  of  near  Tj^  symmetry  is  obtained  regardless  of 
the  /-orbital  interaction. 
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